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ABSTRACT 

We measure the gas-phase oxygen abundances of ~ 3000 star-forming galaxies at z = 0.05 — 0.75 
using optical spectrophotometry from the AGN and Galaxy Evolution Survey (AGES), a spectroscopic 
survey of Jab < 20.45 galaxies over 7.9 deg^ in the NOAO Deep Wide Field Survey (NDWFS) 
Bootes field. We use state-of-the-art techniques to measure the nebular emission lines and stellar 
masses, and explore and quantify several potential sources of systematic error, including the choice 
of mctallicity diagnostic, aperture bias, and contamination from unidentified active galactic nuclei 
(AGN). Combining volume- limited AGES samples in six independent redshift bins and ~ 75, 000 star- 
forming galaxies with tab < 17.6 at z = 0.05 — 0.2 selected from the Sloan Digital Sky Survey (SDSS) 
that we analyze in the identical manner, we measure the evolution of the stellar mass-metallicity 
{Ai — Z) between z = 0.05 and z ~ 0.75. We find that at fixed stellar mass galaxies at z ~ 0.7 have 
just 30% — 60% the metal content of galaxies at the present epoch, where the uncertainty is dominated 
by the strong-line method used to measure the metallicity. Moreover, we find no statistically significant 
evidence that the M. — Z relation evolves in a mass-dependent way for M. ~ 10^ * — 10^"'^ M.q star- 
forming galaxies. Thus, for this range of redshifts and stellar masses the M. — Z relation simply shifts 
toward lower metallicity with increasing redshift without changing its shape. 

Subject headings: galaxies: abundances — galaxies: evolution — galaxies: fundamental parameters 



1. INTRODUCTION 

Like stellar mass, the gas-phase metallicity of a galaxy 
is a sensitive observational diagnostic of its past star for- 
mation history and present-day evolutionary state for 
the simple reason that in a closed system metallicity 
increases monotonically with each successive generation 
of massive starsQ In detail, however, galaxies are not 
closed systems: infall of cold, metal-poor gas from the 
intergalactic medium, rapid gas accretion via minor and 
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^'^ Throughout this paper we use the terms metallicity and abun- 
dance interchangeably to mean the heavy-element content of the 
warm (T 10* K) interstellar medium of galaxies. In particu- 
lar, we make the reasonable assumption that the nebular oxygen 
abundance, historically written as 12 -|- log (O/H), traces the total 
gas- phase metallicity, Zgas- 



major mergers, and supernova-driven winds of gas and 
metals can modulate the metallicity of individual galax- 
ies according to their large-scale environment, gas sup- 
ply, and assembly history. Therefore, accurate abun- 
dance measurements provide valuable insight into the in- 
terplay between many fundamental processes in galaxy 
evolution, including star formation, gas accreti on, and 
supernova-driven feedback acr oss cosmic time (jTinslevI 
ITosnl : [PagelllQQTl : IPettinil[200l . 

Nebular emission lines — ubiquitous in the rest-frame 
optical spectra of star-forming galaxies — provide a par- 
ticularly powerful way to study the chemical abundances 
of both nearby and distant galaxies. Among the most 
commonly observed lines are the Balmer Ha A6563, 
H/3 A4861, and H7 A4340 hydrogen recombination 
lines, and the coUisionally excited [O 11] AA3726, 3729, 
[Ne III] A3869, [O iii] AA4959, 5007, [N 11] AA6548, 6584, 
and [S 11] AA6716, 6731 forbidden lines. Because these 
lines originate principally in star-forming (H 11) regions, 
they trace the physical conditions in the gas from which 
the current generation of massive stars is forming. The 
metal lines in particular are the principal coolants in H 11 
regions, making them sensitive to the total abundance 
of oxygen and other heavy elements in the interstellar 
medium. Moreover, the relative line-strengths can be 
used to infer the interstellar pressure, density, temper- 
ature, ionizing radiation field strength, dust reddening, 
and the presence of an active galactic nucleus (AGN). 
Finally, from an observational standpoint, many of the 
nebular lines are intrinsically strong, making them mea- 
surable in even relatively low signal-to-noise spectra over 
a broad range of redshifts. 

Following the discovery of the correlation between dy- 
namical mass and oxygen abundance in dwarf irregular 
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galax i es (iLeaueux et al.lll979l: s ee also lGarnett fc Shields! 
Il987t lOev fc KennicuttI 119931) . numerous subsequent 
studies demonstrated that star-forming galaxies rang- 
ing from the lowest-luminosity dwarfs to massive 
disk galaxies obey a well-defined luminosity-mctallicity 
(L — Z) correlation spanning several orders of magni- 
tude in i3-band luminosity (iGarnett fc Shields! 0*9871: 
Skillman et"all!l989!: !Zaritskv et al.!!l994 !Pilvugin et ah 



2004|lLamareilie et al.!l2004r!Salzer et all2005!: !Lee et al 
20061: Moustakas et al.ll2010l ). Metalhcity was also found 
to correlate with morphologi cal type (Edmunds Ik Pagcl 
T98l. surface mass density (jMcCall! 11983 iRvdeilTTgga 



Garnett et all!l997!) and max i mum rotational velocit y 
(jZaritskv et al.! !1994!: !Garnett! [200l !Dalcanton! [20071 ). 



However, because these global properties all correlate 
with one another, the underlying physical driver of the 
L — Z relation re m ained elusive. 

!Tremonti etlll (pOOl was the first to leverage 
the tremendous statist ical power of th e Sloan Digi- 
tal Sky Survey fSDSS: lYoFk et al.l!200dl ) spectroscopic 
database to show that the gas-phase metallicity of 
a galaxy correlates best with stellar mass, Ai, now 
known as the stellar mass-mctallicity (Ai — Z) rela- 
tion. The A4 — Z relation reveals that the gas-phase 
metallicity of star-forming galaxies increases monoton- 
ically with stellar mass and then remains relatively 
constant above « 3 x 10^° A^o- Although 

subsequent studies have found weak residual corre- 
lations from the Ai — Z relation with stellar mass 
densitv (!Tremonti et ahl !2004!: !Liang et all !2010t) . size 
Ellison et al.!!2008al). large- and small-scal e environment 



.Mouhcine et al! 12007!: !Cooper et a l]!2008!: lEUi son et al.l 
!2008bl !2009l: !Peeples et al.l !2009!). star formation rate 



(SFR: !Mannucci et al.! 120101: !Lara-L6pez et al! !2010b!: 
! Yates et'al _ 2011 Cresci et al.!!201lD . and the existence 
of bars I Ellison et al.!!201lD . the intrinsic dispersion in 
the Ai — Z relation is < 0.1 dex, making it among the 
tightest empirical correlations known. 

The physical origin of the Ai — Z relation remains 
under debate. One possibility is that low-mass galax- 
ies started forming stars later than massive galax- 
ies (i.e . , th ey are "younger"; !Noeske et al.! !2007a! : 
!LeitneH !2011t ) and have been less efficient at synthe- 
sizing metals via star formation (!Brooks et al.l !2007! : 
!Mouhcine erall!2008!: !Calura et al.l!2009[) ThiThTterpre- 
tation is qualitatively consistent with the measured low 
gas fracti ons in massive galaxies relative to lower-mass 
galaxies (IMcGaueh k de^B \Twh IGeha et al.l l200l 
!Garcia-Appadoo et al.l!2009[) . and with the observed cor- 
relation between stellar mass and SFR ( t he star for- 
mation sequence; 'Brinchmann et ahl 20041: Salim et al.l 
[2007; Elbaz et al. 2007; Nocskc^eLalJ [2007^) . Another 
popular interpretation is that supernova-driven galactic 
win ds preferentia ll y expel metals from low- mass galax- 
ies ([Larson! 
!Dalcantoiil 



19^ iGarnetti \20U% [Tremonti et"aI[l2BP: 
2007[ ). Indeed, state-of-the-art hydrody- 



namic and semianalytic theoretical models require metal- 



erties, including the A4 — Z i 


elation (!Kobavashi et al. 


2007 


; [Finlator & Dave! !2008!: 


Dutton & van den Bosch 


2008 


; !Peeoles & Shankad !2011 


; !Dave et al.l!2011bD. A 



related outstanding question is whether gas and 
metals ejected in previous sta r formation episodes 
are re-accreted onto galaxies ([De Lucia et al.l [2004! : 



!Oppenheimer et"aI1 l2010f ). or whether galaxies are 
predominantly fed by infall of c old, metal-poor gas 
from the intcrgal actic medi u m ([K5ppen fc Edmundi 
19991: [Fauc her-Gimereetall |2011[). Alternatively, 
Koppen et al.l (|20C)7l ) show that a SFR-dependent, and 



therefore mass-dependent, stellar initial mass function 
(IMF) naturally explains the observed A4 — Z relation 
without needing to invoke metal-enriched outflows from 
galaxies. 

One way of gaining insight into the complex and mul- 
tifacetcd interrelationship between chemical abundance 
measurements, star formation, gas accretion, and the 
role of supernova feedback during galaxy growth is to 
measure the evolution of the Ai — Z relation. The time 
since z = 1, spanning ~ 60% of the age of the Universe, 
is important for man y reasons. Measurements of the star 
form ation sequence ([Noeske et all !2007bl: !Wuvts et all 
[201 ll ) suggest that star-forming galaxies over this red- 
shift range evolve smoothly as a population, driven by 
continuous, secular processes like gas consumptio n and in 
situ s tar formation ([Bell et al.ll2005l:lNoeskd2009HLeitiieil 
!2011f ). Bolstered by measurements of their morpholog- 
ical distribution (jBclLet al. 2005.; K pnis hi et al 201J) 
and clustering properties (jCoil et all [20081 : IZehavi et al.l 
l2011f ). these results suggest that we can draw direct 
evolutionary connections linking star-forming galaxies 
across this redshift range. In addition, the metallicity- 
sensitive R23 parameter (defined in ij5|) can be measured 
from ground-based optical spectroscopy of galaxies at 
least to z ^ 0.8, ensuring that metallicities can be de- 
rived using a single consistent abundance calibration. 

Unfortunately, previous studies have reported widely 
varying results on the evolution of both the shape 
and normalization of the Ai — Z relation since z = 
1. [Sayagiio et al. (200^ and [Zahid et all (12011.) find 
that the mean metallicity of galaxies at z '--^ 0.8 with 
M - 10i°-3 A4q differs by just - 0.05 dex (factor of 
~ 1.12) relative to similarly massive star- forming galax- 
ies at z ~ 0.1, whereas they find that galaxies with 
Ai ^ 10^-^ Mq undergo a factor of 2 — 3 more chemical 
evolut ion over the same re dshif t interval. On th e other 
hand, !Liang et'al] (|2006aj) and [Cowie fc Barged ([2008h 
report 0.2 — 0.3 dex (factor of 1.6 — 2) of chemical evolu- 
tion for galaxies with Ai ~ lO^'^ '^ AiQ since z ~ 0.7 and 
no statistically significant evidence for evolution in the 
shape of the Ai — Z re l ation (see also iRodrigues et al 



[2008D . !Lamareille et al.! (|2009! ) and !Perez-Montero et al 
(|2009( ) find an even greater amount of metallicity evolu- 
tion, ^ 0.35 dex (factor of 2.2) since z = 0.6 — 1, and 
a flattening of the Ai — Z relation at intermediate red- 
shift. Other studies based on the evolution of the B- 
band L — Z relation report that at fixed luminosity star- 
forming galaxies at z = 0.5 — 1 are 0.1 — 0.7 dex (fac- 
tor of 1.2 5 — 5) more metal -poor t han l o cal sta r-forming 
galax i es ([Kobulnickv & Z aritskv! !1999!: !Caroll o & LilH 
200I ILiUv et all 120031: iKobulnickv fc Phillips! pOOlF 
Kobulnickv fc KewlevI 120041: IMaier et all 12004 HooT 



2006t iLamareille et al.l 12009! ). However, the signifi- 



cant amount of luminosity evolution e xperienced by 
blue, star-forming galaxies since z = 1 ([Blanton! [2001 : 
[Faber et all!2007! : (Cool et al.!!201l! ) renders the interpre- 
tation of the L — Z relation less straightforward. 
The origin of these widely varying results on the cvo- 
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lution of the M. — Z relation can be attributed to vari- 
ous issues, including cosmic variance, small sample size, 
heterogenous selection criteria, systematic differences in 
the methods used to infer nebular abundances, an in- 
consistent analysis of the local M — Z relation, spec- 
troscopy with insufficient spectral coverage, signal-to- 
noise (S/N) ratio, or instrumental resolution, and uncer- 
tain stellar mass estimates due to limited broadband pho- 
tometric coverage. Consequently, the efficacy of existing 
observational metallicity constraints on state-of- the-art 
theoretical models of galaxy formation (jPe Luci a ^ alJ 
200llBrooks et al.ll2007l : iDutton fc van den Boschll2009l : 



Dave et al.ll2011bl )has been fairly limited. From an ob- 



servational standpoint, there exists a clear need for a 
large, homogeneously selected sample of galaxies with 
high-quality spectroscopy and reliable stellar mass esti- 
mates to better constrain the evolution of the M. — Z 
relation at intermediate redshift. 

To address this need, we measure the evolution of the 
M. — Z and 5-band L — Z relations at intermediate red- 
shift using oxygen abundances of ~ 3000 star-forming 
galaxies at z = 0.05 — 0.75 observed as part of th e AGN 
and Galaxy Evolution Survey (AGES:lKochane k et al.l 
120 lit ), and - 75, 000 galaxies at z = 0.05 - 0.2 selected 
from the SDSS. The AGES main galaxy survey con- 
sists of optical spectrophotometry for ~ 12,000 galax- 
ies in the ~ 9 deg^ NOAO Deep Wide F ield Survey 
(NDW FS: ITannuzi fc Dex^llOQl iBro^ et al.l [2003. 2007, 
|2008| ) Bootes field at a median redshift of z ~ 0.3. This 
sample is statistically c omplete over 7.9 deg^ for galax- 
ies with /ab < 20.45 (jKochanek et al.l[2011[ ). reaching 
^ 2 mag deeper t han the SDSS main galaxy sample 
(jStrauss et al.ll2002| ) over a considerably larger area than 
other recent or ongoing sur veys of intermediat e-redshift 
galaxies su ch as DEEP2 (iDavis et al.l 120031) . AEGIS 
(|Davis et a l. 2007), VVDS (|Le Fevre et"aLn20 04, 2001; 
IGarilh et al..,2008.1 . and zGOSMOS (ILillv et al..,2009.l FI 
The statistical completeness, large sample size, and avail- 
ability of high-quality optical spectroscopy and deep op- 
tical and near-infrared photometry enables us to investi- 
gate the evolution of both the M — Z and optical L — Z 
relations of star-forming galaxies at z = 0.05 — 0.75 us- 
ing a single, self-consistent abundance diagnostic over the 
entire redshift range. 

The plan of the paper is as follows. In Sj2]we present the 
ground-based optical and near-infrared imaging we use, 
summarize the AGES observations, and describe how we 
measure the nebular emission lines for the galaxies in our 
sample. In fJ3]we describe the methodology used to derive 
rest-frame luminosities, colors, and stellar masses, and in 
21 we select a subset of the AGES galaxies for chemical 
abundance analysis. We present the methods we use to 
derive oxygen abundances in ^ and our principal results 
in fjni where we quantify the mass-dependent evolution of 
the M — Z relation for star-forming galaxies since z = 
0.75. In ^ we address the effect of various potential 
sources of systematic uncertainty on our results, and in 
^ we discuss recent theoretical work on the origin and 
evolution of the M. ~ Z relation. Finally, we summarize 

For galaxies at z < 0.5, the Galaxy and Mass Assembly 
(GAMA) survey will ultimately supplant all these surveys by ob- 
taining optical spectroscopy for several h undred thousand g alaxies 
brighter than tab ~ 19-8 over 360 deg^ HDriver et al.ll20lffl . 



our principal conclusions in ^ 

We adopt a concordance cosmology with 51,„ ~ 0.3, 
rjA = 0. 7, and feyn = ffn/l OO = 0.7 the AB magni tude 
system (|Oke fc Gunnlll983[ ). and the iChabrieJ (|2003l ) ini- 
tial mass function (IMF) from 0.1 — 100 M.q unless oth- 
erwise indicated. For reference, the conversion from Vega 
to AB for the /-band filter used to select AGES targets 
is -1-0.45 mag, that is /ab = /voga + 0.45. 

2. OBSERVATIONS 

In the following sections we present the multiwave- 
length imaging of the Bootes field that we use ( §2.ip . 
describe the AGES optical spectroscopy and emission- 
line measurements (i i2.2p . and construct a comparison 
sample of local galaxies from the SDSS {\ 



2.1. Multiwavelength Photometry 
2.1.1. Optical and Near-mfrared Imaging 

Our baseline optical observations consist of deep 
BwRI imaging available as part of the NDWFS third 
public data releaseE3 The imaging was carried out us- 
ing the MOSAIC-I camera at the KPNO/MayaU 4 m 
telescope, reaching a 5cr depth of ^ 26.5, ^ 25.5, and 

25.3 mag in a 2" diameter circular aperture in /3vVi 
and /, respectively. Note that although A"-band imag- 
ing from the NDWFS is also available for ^ 60% of the 
Bootes survey region, we do not use these older data in 
favor of the more recent near-infrared (near-IR) observa- 
tions described below. 

We supplement the NDWFS photometry with U - and 
z-band observations obtained as part of two other imag- 
ing campaigns. F. Bian et al. (2012, in prep.) have 
obtained deep J7-band imaging of the Bootes field us- 
ing t he Large Binocular Gamera (LBG; iGiallongo et al.l 
I2008D mounted at the primary focus of the 8.4 m Large 
Binocular Telescope (LBT); these observations reach a 
5(7 depth of ^ 25. 2 mag for p oint sources. Second, the 
zBootcfH survey (|Cooll 120071 ) imaged 7.62 deg^ of the 
Bootes field to a 3(t depth of 22.4 mag in a 3" diameter 
aperture using the QOprime prime-focu s wide-field im- 
ager mounted at the Bok 2.3 m telescope (jWilliams et al.l 
[200l . 

Finally, we utilize deep JHKg near-IR imaging of the 
entire Bootes field obtained using the NOA O Extremely 
Wide -Field Infrared Mosaic (NEWFIRM; lAutrv et al.l 
|2003[ ) instrument at the KPNO/MayaU 4 m telescope. 
The NEWFIRM observations achieved a 5a depth of 
~ 22.9, ^ 22.1, and ^ 21.4 mag in a 3" diameter aper- 
ture in J, H, and A'^, respectively. More details regard- 
ing the NEWFIRM observations and reductions will be 
described in an upcoming paper (A. Gonzalez, 2011, pri- 
vate communication). 

2.1.2. Ancillary X-ray, Mid-Infrared, and Radio Imaging 

We use observations of the Bootes field in the X-ray 
from the Chandra X-ray Obs ervatory, at 3.6 — 8 //m from 
the Spitzer Space Telescope (| Werner et al.|[200^ . and at 
1.4 GHz from the Westerbork Synthesis Radio telescope 

The BiY filter is a "wide" filter with a bluer effective wave- 
length (A£^^_2_4200 A) than the standard Johnson-Morgan B-band 
filter IIBesselllll99^ 

http : //archive .noao ■ edu/nsa/zbootes .htmll 
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to help identify and remove AGN from our sample (see 
^4.2|) . Here, we briefly describe these ancillary data. 

The 5 ks X-ray observations were obtained with 
the Chandra Advanced CCD Imaging Spectrometer 
ACIS) instrument as part of the XB ooteJ*^ survey 
Murrav et alll2005bl : iKenter et alJl2005( ). XBootcs cov- 
ered ~ 8.5 deg^ of the Bootes field to an on-axis lim- 
iting flux of ~ 7.8 X 10~^^ erg cm~^ in the 
0.5 — 7 keV band. We use the publically available 
catalog of 3213 X-ray point-sources with four or more 
X- ray counts matched to the NDWFS optical catalog 
bv lBrand etHI (pOOl using a Baycsian matching algo- 
rithm. 

In the i nfrared, the Spitz er Deep Wide-Field Survey 
(SDWFS: lAshbv et al.l l2009t) obtained deep, muhi-epoch 
mid-IR imaging of the Bootes field with the Sp itzer In- 
frared Array Camera (IRAC; iFazio et al.ll2004D . These 
observations reach aperture-corrected 5a depths of 
22.6, - 22.1, - 20.2, and ~ 20.2 mag in a 4" diameter 
aperture at 3.6, 4.5, 5.8, and 8 /im, respectively. In our 
analysis we use the 3.6 /im-detcctcd, aperture-corrected 
4" diameter 3.6 — 8 /im aperture magnitudes available 
at the NA SA/IPAC Iri f rarcd Science ArchiveO, as de- 
scribed bv lAshbv etHI (|2009( ). 

Finally, radio observations from the WSRT 1.4 GHz 
radio survey covered ~ 7 d eg^ of the Bootes fi eld to a 
5(7 limiting flux of 140 fxJy (|de Vries et al.ll2002f ). When 
matching to t he AGES catalog we use a 3" diameter 
search radius (jHickox et al.ll2009| ). 

2.2. AGES Optical Spectroscopy 

The AGES main galaxy survey targeted galaxies 
brighter than Iab = 20.45 and covered 7.9 deg^ of 
the NDWFS Bootes field, making it among the largest 
wide-area spectroscopic surveys of intermediate-redshift 
galaxies conducted to date. We briefly summarize the 
AGES observations in t|2.2.11 and in ^2.2.2\ we describe 
the procedure used to measure the nebular emission-line 
strengths from these data. A more thorough descrip- 
tion of the AGES experimental design, o bservations, and 
redshift m easurements can b e found in iKochanek et al.l 
((20T1I) and lCool et all (IMl . 

2.2.1. Observations, Reductions, and Spectroscopic 
Completeness 

The MMT/Hcctospec multi- fiber optical spectrograph 
is fed by 300 robotically controlled 1'.'5 diameter fibers, 
enabling efficien t data acquisition over a 1° diame- 
ter fi eld-of-view (jRoh et al.lll998t iFabricant et al.lll998l 
[20051) . The 270 line mm"! grating used by AGES 
provides 3700 — 9200 A spectra with 1.2 A pixels at 
~ 6 A FWHM resolution. An atmospheric dispersion 
corrector (ADC) built into the MMT f/5 wide-field lens 
eliminates wavelength-dependent light-loss due to atmo- 
spheric refraction, while relative spectrophotometric cal- 
ibration is facilitated by observin g F sub-dwarf standard 
stars selected from the SDSS (jAbazaiian et al.l 120041 : 
IFabricant eralll2008[) . 

We reduce the AGES spectra using standard pro- 
cedures implemented in HSREeF^. a customized Hec- 

^® http://wMW.noao . edu/noao/noaodeep/XBootesPublic 
1'' http://irsa.ipac.caltech.edu/data/SPITZER/SDWFS 
http : //www . astro . princeton . edu/$\sim$rcool/hsred] 



tospec data reduction package based on the SDSS 
spectroscopic data reduction pipelincEEl Briefly, we 
bias- and overscan-subtract the data, and then extract, 
wavelength-calibrate, flat-field, and sky-subtract each 
galaxy spectrum. We combine multiple exposures us- 
ing invers e variance weights, and robustly reject cos- 
mic rays (|van Dokkum|[200lr) . We then flux-calibrate 
each spectrum, divide by a suitably scaled telluric ab- 
sorption spectrum, and correct for foregro und Galac- 
tic reddening (Ry = Avj E(B - V) = 3.1; lO'Donnelll 
IT99I ISchlegel et aD[T998h . Note that five of the Hec- 
tospec configuration obtained in 2004 could not be 
flux-calibrated, and have been excluded from the present 
analysis. Finally, we determine redshifts using two in- 
dependent cod es and cross- validate t he results by visual 
inspection (see IKochanek et al.ll201l] for details). 

Because we are interested in faint emission lines at in- 
termediate redshift, accurate sky subtraction is crucial. 
Therefore we implemented an automated principal com- 
ponent analysis technique originally developed for the 
SDSS to suppress the amplitude of the sky-subtraction 
residuals redward of ^ 6700 A, which a re dominated by 
a "fo rest" of atmospheric OH sky lines (jWild fc Hewettl 
l2005f l. Briefly, the method reconstructs a model of the 
night-sky spectrum for each configuration of the Hec- 
tospec fibers by exploiting the fact that the OH sky lines 
are highly correlated both in wavelength and between all 
the fibers. On average, we obtain a factor of two im- 
provement in the S/N ratio of affected pixels. 

The AGES spectra suffer from two additional data- 
reduction issues that do not affect the redshift measure- 
ments, but which are relevant to our emission-line analy- 
sis. First, the ADC was operated incorrectly during some 
of the observations that took place in 2004, leading to 
systematic errors in the spectrophotometry blueward of 
^ 50 00 A for these early observations (jFabricant et al.l 
l2008f ). And second, the spectra are occasionally con- 
taminated by the featureless continuum of one or both 
of the light-emitting diodes (LEDs) located on the Hec- 
tospec fiber-positioners. In most cases this extra light 
appears as a rapidly rising continuum redward of 8500 A, 
although it can also bias the average sky spectrum, 
resulting in an oversubtraction of the sky fiux in the 
8500 — 9200 A wavelength range. We address these issues 
in g222I 

We conclude this section with a brief discussion of 
the AGES spectrosco pic completeness (fo r details see 
IKochanek et aIll20Tll and lCool et al.ll201lh . AGES cov- 
ered 7.9 dcg^ of the Bootes field with 15 slightly overlap- 
ping pointings, each with three configurations. Within 
each field, AGES targeted 100% of galaxies with Jab < 
18.95 and employed a variety of sparse-sampling criteria 
to observe a subset of galaxies with 18.95 < /ab < 20.45 
that were also bright at other wavelengths. Because these 
sampling fractions are known exactly, it is easy to recover 
a sample that is statistically complete to /ab < 20.45 by 
simply weighting by the inverse of the sampling rate. 
The other sources of incompleteness in AGES are rela- 
tively small. Approximately 5% of objects in the parent 

|http : //spectro ■ princeton ■ edu] 

A configuration refers to a given geometrical layout of the 300 
Hectospec fibers. 
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sample failed to be assigned a spectroscopic fiber {fiber 
incompleteness) . Among the objects that were observed, 
AGES failed to measure a redshift for ^ 2% {redshift in- 
completeness) , which is a weak function of surface bright- 
ness. Finally, we estimate that roughly 4% of objects are 
missing from the targeting catalog [photometric incom- 
pleteness) due to proximity to a bright star or problems 
with the photometry itself. To correct for these effects, 
we weight every object in AGES by the product of these 
terms, resulting in a sample that is statistically complete 
to /ab = 20.45. 

2.2.2. Emission-line Measurements 

High-resolution population synthesis models have 
become an indispensible tool for investigating the inte- 
grated optical spectra of g alaxies fBruzual & Chariot 
200l ivkzauez fc Leithereif [200 5: GriUo et al. 200, 



Conrov et al.ll2009t iPercival et al. 2009: Vazdekis et al 
2010[) . Bv modeling and subtracting the stellar contin 



uum from the observed spectrum, these models have 
made it possible to study the optical emission lines free 
from the systematic effects of Balmcr and metal-line ab- 



sorption (e.g., iP anter et al 
Cid Fernandes et all 120051: 



2006 




2007 


: 



2003; Trcmonti et al. 20a 
]]\'Ioustakas fc Kennicut 



Sarzi et al.ll2006l: lOcvirket al.,.200& : ,Toieiro et al 
Oh et al.ll2011[ ). Neglecting the effects of stellar 



absorption, particularly under the Balmer emission lines, 
can lead to severe biases in the nebular abundances, 
SFRs, and dust re ddenings inferred from the optical 
emission lines (e.g.. iKennicuttI 119921: iKobulnicky et al.l 

■ m 



1999; 



201C; 



Rosa-Gonzaiez et al.ll2002t iMoustakas et all 12001 
Asari et al.ll2007[) . 



Our basic strategy is to model the observed spectrum 
of each galaxy as a non-negative linear combination of 
simple (i.e., instantaneous-burst) population synthesis 
models of varying ages, making some simplifying assump- 
tions regarding the stellar metallicity of the galaxy and 
the effects of dust attenuation. We then subtract the 
best-fitting continuum model from the data, and fit the 
residual emission-line spectrum assuming Gaussian line- 
profiles. We emphasize that our principal goal is to ob- 
tain an emission-line spectrum that has been corrected 
self-consistently for stellar absorption, not to constrain 
the star formation and chemical evolution history of the 
galaxy from its integrated stellar spectrum. The multi- 
dimensional parameter space of star formation history, 
age, metallicity, and dust attenuation is highly degener- 
ate (i.e., there are many local minima), and our simple 
ap proach is ill-suited fo r finding the global solution (but 
see lWalcher et aLll2010l and references therein, for more 
advanced techniques that tackle this issue) . Fortunately, 
the emission-line strengths we measure do not depend 
sensitively on the simplifying assumptions underlying our 
technique. 

We construct our template set using the 
iBruzual fc CharlotI ([200l hereafter BC03) popula- 
tion synthesis models, based on the Padov a 1994 stellar 
evolutionary tracks (jGirardi et al.l I1996L and refer- 
ences therein) a nd the empirical STELIB stellar library 
(|Le Borgne et a l. 2003). We choose 10 instantaneous- 
burst, solar-metallicity [Z = 0.02) models with ages 
spaced quasi-logarith mically in time b etween 5 Myr and 
13 Gyr, assuming the lChabrieJ ()2003f ) initial mass func- 
tion (IMF) from 0.1 — 100 Mq. We treat dust reddening 



as a free parameter, and adopt the iCalzetti et all (|20000 
dust attenuation law. Finally, we find the best-fitting, 
non-negative linear combination of the templates using a 
modified version of the pPXII^ continuum-fitting co de 
(|Cappellari fc Emse Ucnr 2004: M oustakas et al.ll2010( ). 

We verified that adopting different model parameters 
does not have a significant effect on the derived emission- 
line strengths. Specifically, we experimented with a 
wider range of stellar metallicities, varied the number of 
templates (i.e., instantaneous-burst ages), tried several 
different extinction curves, and assumed a different ini- 
tial mass function. These results indicate that it is not 
necessary to find the global solution to the continuum 
fitting problem in order to be able to study the emission- 
line properties of galaxies free from the systematic effects 
of stellar absorption. 

We fit the observed stellar continuum of each galaxy it- 
eratively. First, we isolate the rest-frame ^ 3600—4400 A 
spectral range to derive a more precise measurement of 
the absorption-line redshift, Zabs, and to obtain an esti- 
mate of the stellar velocity dispersion, Cdisp, accounting 
for the instrumental resolution of our spectra (^ 6 A 
FWHM) and the resolution of the BC03 models (~ 3 A 
FWHM; BC03). Next, we fix Zabs and a^isp at the de- 
rived values and model the full wavelength range, aggres- 
sively masking pixels that might be affected by emission 
lines, sky-subtraction residuals, telluric absorption, or 
the red leak (see §2.2.ip . For approximately one-quarter 
of the observed Hcctospec configurations (largely from 
2004; see ij2.2.ip there are non-negligible errors in the 
spectrophotometry shortward of Aobs ~ 5000 A. There- 
fore, for these configurations, we refit the spectra, mask- 
ing out the affected blue observed- frame wavelengths, 
and then divide the observed spectrum by its best-fitting 
model, in effect constructing a sensitivity function cor- 
rection for each galaxy. We then compute the median 
sensitivity function correction using all the spectra on 
a given configuration, and divide all the spectra by the 
median correction. Finally, we refit those galaxy spectra 
one last time, this time using the full wavelength range. 
We assess the effect of this correction on our abundance 
analysis below. 

The last step before fitting the emission lines is to sub- 
tract the red leak and any remaining low-level residuals 
due to imperfect sky-subtraction, template mismatch, or 
imperfect flux-calibration. Given the best-fitting contin- 
uum model for each galaxy, we subtract it from the data 
and remove these residual differences (typically of order 
a few percent) using a sliding 151-pixel median filter. 
Finally, we subtract this median residual spectrum from 
the original data, and then refit and subtract the stel- 
lar continuum one last time, leaving a pure emission-line 
spectrum. 

Following [Moustakas et all (|2010D . we fit the emission- 
line spectrum of each galaxy iteratively using a 
modified version of the GANDALi|3 emiss ion-line fit- 
ting c ode, assuming G aussia i i line-profiles ([Sarzi et al.l 
[20M : ISchawinski etaD [20071: lOh et al.l I201ll) . Soecif- 
ically, we fit the [Ne v] A3426, [O ii] A3727, 
[Ne III] A3869, [O iii] AA4959,5007, [N ii] AA6548,6584, 

'http : //www- astro .physics . ox. a c .uk/SNsimSmxc/ldl | 
htt p : //www . strw. leidenun iv. iil/sauron| 
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Fig. 1. — Comparison of EWs measured from the fluxed and 
unfluxed AGES spectra for H/3, [O III] A5007, and [O ll] A3727. 
In each panel the solid line shows the one-to-one relation, and 
the contours enclose 50% and 95% of the points. We find ex- 
cellent statistical agreement among the two emission-line mea- 
surements, assuming that EW(H/3) measured from our unfluxed 
spectra is subject to 2.8 ± 1.5 A of stellar absorption (solid and 
dashed red lines in the upper-left panel). In the lower-right panel 
we compare the metallicity-sensitive EW(-R23) = (EW([0 ll]) -f 
EW([0 lll]))/EW(H/3) parameters measured from the fluxed and 
unfluxed spectra, after statistically correcting EW(H/9) measured 
from the unfluxed spectra for an average 2.8 A of stellar absorption. 
These comparisons indicate that in spite of the spectrophotometric 
errors in AGES, our oxygen abundance estimates are not adversely 
affected. 

and [S ii] AA6716, 6731 forbidden lines, and the first 
five lines in the Balmer series: Ha A6563, H/3 A4861, 
H7 A4340, H(5 A4101, and H5 A3889. On the first iter- 
ation we tie the redshifts and intrinsic velocity widths 
(accounting for the instrumental resolution) of all the 
lines together to aid in the detection and deblending of 
weak lines. In addition, to reduce the number of free 
parameters, we constrain the [O ill] A5007/[O ill] A4959 
and [N 11] A6584/[N 11] A6548 doublet ratios to be 3:1 
(jOsterbrock fc Ferla nd 2006). On the second iteration 
we relax most of these constraints and use the best-fitting 
parameters from the first iteration as initial guesses. This 
second step is necessary because of our imperfect knowl- 
edge of the wavelength-dependent instrumental resolu- 
tion, and to account for the fact that [O 11] A3727 is a 
doublet which is better represented at the spectral reso- 
lution of our data as a single, slightly broader Gaussian 
line than two closely-spaced Gaussian line-profiles. Note 
that even on the second iteration, however, we (sepa- 
rately) constrain the redshifts and velocity widths of the 
[O III], [N 11], and [S 11] doublets to have the same values. 

We evaluate the significance of each emission line by 
comparing the amplitude of the best-fitting Gaussian 
model, Al , against the standard deviation of the residual 
(continuum-subtracted) spectrum, CTc, around the line; 
we define a significant line as having Ai > 3ac- For each 
undetected emission line we estimate the la upper limit 



on the total line-flux assuming a Gaussian line-profile 
using 

F = V2^a,S,, (1) 

where Sy is the mean velocity width of all the detected 
emission lines. Finally, we estimate the equivalent width 
(EW) of each emission line by dividing the integrated 
flux of the Gaussian model by the mean flux of the stel- 
lar continuum immediately blueward and redward of the 
given line. 

We conclude this section by evaluating the effect of 
spectrophotometric errors on our line-measurements by 
measuring the [O 11], H/3, and [O ill] A5007 emission- 
line EWs from the unfluxed AGES spectra. For these 
spectra we model the stellar continuum as a low-order 
B-spline with iterative outlier rejection, subtract it from 
the data, and measure the EW of each emission line as 
described above. In Figured] we compare the [O 11], H/3, 
and [O III] emission-line EWs measured from the fluxed 
and unfluxed spectra, as well as the metallicity-sensitive 
EW(i?23) = (EW([0 ii])-f-EW([0 iii]))/EW(H/3) pa- 
rameter (see §5.2p . We render the data in greyscale with 
the contours enclosing 50% and 95% of the points, and 
plot the one-to-one relation as a solid line. Overall we 
find good statistical agreement among the EWs mea- 
sured from the two sets of spectra, although H/3 is clearly 
affected by stellar absorption. However, applying a mean 
2.8 ± 1.5 A correction to the EW(H/3) values measured 
from the unfluxed spectra, shown as solid and dashed red 
lines in the upper- left panel of Figure [1] would bring the 
two EW measurements into good agreement. Meanwhile, 
the statistical agreement for EW([0 11]) and EW([0 ill]) 
is excellent. In the lower-right panel we compare the two 
EW(i?23) estimates after correcting EW(H/3) measured 
from the unfluxed spectra for an average 2.8 A of stellar 
absorption. These comparisons indicate that our oxygen 
abundance estimates are not adversely affected by the 
spectrophotometric errors in AGES. 

2.3. SDSS Comparison Sample 

At low redshift (z ~ 0.1) AGES is signiflcantly affected 
by variations in large-scale structure (i.e., cosmic vari- 
ance) because of the relatively small cosmological vol- 
ume probed. Therefore, we use observations obtained 
as part of the SDSS as a local benchmark for studying 
the chemical abundance properties of the more distant 
galaxies observed by AGES. 

The SDSS has obtained ugriz broadband photometry 
and high-precision optical (3800 — 9200 A) spectropho- 
tometry for roughly t wo million object s over nearly one- 
quart er of the sky ()York et al.l 120001 : iStoughton et al.l 
[2OOI FI We draw our SDSS comparison sample 
from t he NYU Value-Adde d Galaxy Catalog (NYU- 
VAGC; IBlanton et a l.l l2005D corresponding to the SDS S 
Data Release 7 fDR7: lAdelman-McCarthv et al.|[2008[) . 
From this database we select 427, 504 galaxies over 
6955 deg ^ that satisfy the main sample criteria de- 
flned by iStrauss et al.l (I2002D. and l have Galactic ex- 
tinction corrected (jSchlegel et al.lll998D Petrosian magni- 
tudes 14.5 < r < 17.6, and redshifts 0.05 < z < 0.2. We 

http : //www. sdss . org| 

^* htt p : //sdss .ph ysics .nyu. edu/vagc I 
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choose the lower redshift limit to ensure that [O ii] A3727 
lies w ithin the SDSS spectral range (jStoughton et alj 
I2002D . while the upper redshift limit rejects < 2% of 
galaxies in the SDSS main sample. From the VAGC 
we retrieve for each object the spectroscopic redshift, 
the Galactic extinction- corrected uqriz Petro sian and 
cmodel (i.e., total; see lAbazaiian et al.l l2004f ) magni- 
tudes, integrated JHKs photometry (available for ~ 
55% of the sample) from the Two Micron A ll Sky Survey 
Extended Source Catalog (2MASS/XSC; iJarrett et at] 
HOOO; Skrutskie et al. 2006), and the statistical weight 
for each galaxy that corrects for the ^ 8% spectro scopic 
incompleteness of the survey (iBlanton et al.ll2005h . Cor- 
rected for incompleteness the effective number of galaxies 
in our SDSS sample is 464, 145. 

We cross-match this sample against the public 
database of spectroscopic measurements for SDSS/DR7 
galaxies available at the MPA-JHU team websiteI3 The 
MPA-JHU database provides fluxes and EWs for all the 
optical emission lines of interest after carefully model- 
ing and subtracting the stellar continuum of each galaxy 
using an updated, but as-yet unpu blished version of the 
BC03 population synthes i s models (iTreinonti et al.|[2004l : 
iBrinchmann et al.l |20P: iBruzuall 'OOOW . A" negligible 
number of objects, roughly 0.3% of galaxies in the parent 
VAGC sample, do not appear in the MPA-JHU database, 
and have been excluded. We post-process this database 
by adding a flag indicating the presence of a significant 
emission line based on the amplitude of the line relative 
to the variance of the continuum as described in i)2.2.2[ 
and by computing upper limits for each line using equa- 
tion dl]). 

As a consistency check, we retrieved the SDSS spectra 
of a randomly selected sample of 500 galaxies from the 
SDSS data archivj^ and modeled them using the contin- 
uum and e missio n-line fitting code used to fit the AGES 
spectra in §2.2.21 In general the measured line-fiuxes and 
EWs are consistent within the statistical uncertainties. 

3. SPECTRAL ENERGY DISTRIBUTION 
MODELING 

Accurate rest-frame colors, luminosities, and stellar 
masses require photometry that has been measured over 
the same physical aperture for each galaxy. In ^3. II we de- 
scribe the procedure we use to build aperture-matched 
optical to near-IR SEDs of the galaxies in our sample, 
and in ^3.21 and ^ 33.31 we use these data in conjunction 
with the spectroscopic redshifts from AGES to derive 
A'-corrections and stellar masses for the full sample, re- 
spectively. 

3.1. Aperture- Matched Photometry 

Our general strategy is to determine the apparent col- 
ors of each galaxy using aperture photometry measured 
from PSF-matched images, and to use an estimate of the 
total magnitude in the /-band to set the overall normal- 
ization of the SED. This procedure is appealing because 
aperture colors in general have a higher S/N ratio than 
the total magnitude in each band, and it minimize con- 
tamination from neighboring sources. 

^5 |http : / /www ■ mpa-garch lng ■ mpg . de/SDSS/DR7 1 
http : //www . sdss . or g/dr7 | 



To detect sources and to obtain an estimate of the to- 
tal J- band magnitude of each galaxy we run SEXTRAC- 
TOR (iBertin fc Arnoutslll996D in single- image mode on 
the unsmoothed /-band mosaics. We use the mag_AUTO 
(Kron-like) magnitude, /autOi as an estim ate of the in- 
tegrated flux of each galaxy (jKronI ll98Gf ) . Note that 
although MAG_AUTO may miss a significant amount 
of flux at faint apparent magnitudes, and for galax- 
ies with very extended surface-brightness profiles (e.g., 
I Graham fc Driveil l2005| ) , by inserting artifical galaxies 
into the N DWFS /-band mos aics and rc-running S Ex- 
tractor, iBrown et al.l ()2007l ) find that /auto recovers 
the total (input) galaxy magnitude to within ^ 5% for 
galaxies brighter than /ab ~ 21, which is sufficient for 
our purposes. 

Unfortunately, /auto is occasionally biased toward 
bright magnitudes by the low surface-brightness tails of 
bright stars. Although A GES intentionally excl uded re- 
gions around bright stars (|Kochanek et al.|[2011[ ). /auto 
appears to be affected for roughly 10% of the sample. Be- 
cause these low surface-brightness halos are not present 
in the /?-band mosaics, we re- run S EXTRACTOR on the 
/?-band images and use the /?-band mag_AUTO magni- 
tude to obtain a second estimate of the total /-band 
magnitude, Ir^ using the R — I color measured from the 
PSF-matched images in a 6" diameter aperture (see be- 
low). These aperture colors are considerably less affected 
by this ext e nded, low surface-brightness light. Following 
ICool et al.l(l20Tll) . we obtain a new estimate of the total 
/-band magnitude, /', using a statistic that chooses the 
fainter of In and /auto if they differ significantly, and 
averages them otherwise. For ~ 7% (~ 3%) of the sam- 
ple the magnitude correction computed this way exceeds 
0.1 mag (0.5 mag). However, for the remainder of the 
sample there is a very tight correlation between /' and 
-^AUTO: with a median difference of < 0.01 mag and a 
1(7 scatter of ~ 2%. To account for this small amount of 
scatter, in will restrict our statistical (flux-limited) 

galaxy sample t o /jp < 20.4, wh ich excludes ~ 2% of the 
parent sample (jCool et al.ll2011| ). Hereafter, we adopt /' 
as the total /-band magnitude of each galaxy and for 
clarity drop the prime superscript. 

Next, we smooth the optical and near-IR mosaics to a 
common spatial resolution to ensure that the apparent 
colors a re measured from t he same physical area on each 
galaxy (| Brown eTaDIlQOl). The mosaics with the best 
spatial resolution (seeing) are the Bw-, /-, H- and 
A's-band mosaics, so w e smoot h these such that the re- 
sulting stellar PSF is a iMoffatl (|1969f ) profile with 1'.'35 
FWHM and /? = 2.5. The seeing of the U-, z-, and J- 
band mosaics is slightly worse, so we convolve these mo- 
saics to an equivalent 1'.'6 FWHM Moffat PSF. We then 
measure aperture photometry in fixed apertures rang- 
ing in diameter from 1" to 20" centered on the /-band 
position of each source using customized software. We 
use the SExtractor segmentation maps to account for 
pixels contaminated by neighboring objects and missing 
pixels, and estimate the photometric error in each aper- 
ture by computing the interval that encompasses 68% 
of the fiux measured in 100 blank-sky apertures placed 
randomly within 2' of each source. Finally, we construct 
the SED of each galaxy using the 4" diameter aperture 
magnitudes, and then scale everything by the difference 
between the /-band aperture magnitude and the total /- 
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band magnitude. We use 4" aperture magnitudes as a 
compromise between choosing the largest possible aper- 
ture to maximize the S/N, while simultaneously mini- 
mizing contamination from neighboring objects. How- 
ever, we did verify that using the 6" diameter aperture 
magnitudes yields the same colors to within < ±4% and 
with no statistically significant systematic differences. 

3.2. K -corrections 

To estimate the rest-frame colors and luminosities of 
the galaxies in our sample, we compute A'-corrections 
using the publical l y ava ilable code kcorrect^ (v4.1.4; 
iBlanton &: Roweii |2007[) . For each galaxy kcorrect 
takes as input the redshift, the observed photometry, and 
the corresponding filter response curves, and fits the data 
with a non-negative linear combination of five basis tem- 
plates representing a diversity of galaxy star formation 
histories. We have verified that the more detailed SED 
modeling performed in ij3.3l yields the same rest-frame 
quantities, but adopt kcorrect here because of its ease- 
of-use and speed. 

For our AGES sample we fit the aperture-matched 
UBwRIzJHKs photometry assembled in §3.11 and for 
our SDSS sample we fit the observed ugrizJHKg pho- 
tometry (see ^2.3p . all corrected for fore ground Galac- 
tic ex tinction {Ry = 3.1: lO'DonneIilll994l: ISchlegel et al.l 
Il998| ). We compute AT-corrections for each galaxy on a 
case-by-case basis by selecting the bandpass that mini- 
mizes the tr ansformation from the observed frame to the 
rest frame (|Hogg et al.l 120021: IBlanton et al.l [20031 ). For 
example, for our AGES sample we compute the abso- 
lute B-band magnitude (lBesselllll990D from the observed 
Byv magnitude at z < 0.2, from the i?-band magnitude 
at 0.2 z < 0.65, and from the /-band magnitude at 
z > 0.65. We use the best-available filter curve for each 
bandpass by taking into account the detector quantum 
efficiency, telescope throughput, and atmospheric trans- 
mission. Before fitting, we add a minimum photometric 
uncertainty of 0.02 mag in BwRIHKs in quadrature to 
the statistical error in each bandpass. To account for 
the slightly larger PSF of the z-, and J-band mo- 
saics (see t^3.1|) . we assume a minimum error of 0.05 mag 
in these bands. For the SDSS we adopt minimum er- 
ro rs of 0.02 mag in q ri, 0.05 mag in u, and 0.03 mag in 
z (jlvezic et al.l 120041 ). Finally, we also adjust the SDSS 
ugriz magnitudes by -0.036, 0.012, 0.010, 0.028, and 
0.040 mag, respectivel y, to place them on an absolute 
AB magnitude system (jBlanton fc Roweiil2007[) . 

After fitting our AGES sample, we identified small 
residual, systematic ofi^sets between the observed pho- 
tometry and the ma gnitudes syiithesiz ed from the best- 
fitting models (e.g.. Illbert et al]l2005| ). These relative 
zeropoint differences can arise from errors in the pho- 
tometric calibration, slightly mismatched physical aper- 
tures (e.g., due to imperfect PSF matching), imperfect 
filter response curves, and template mismatch. There- 
fore, we fit the sample iteratively, each time adjusting 
the observed photometry in each bandpass by the me- 
dian zeropoint difference between the measured and syn- 
thesized magnitudes. The zeropoints converged to better 
than 0.5% in five iterations. The resulting zeropoint cor- 
rections we apply to the data, all relative to the /-band, 

^'^ |http: //howdy .physics .nyu. edu/ index .php/Kcorrect | 



are: -0.073, -0.032, -0.016, -0.037, -0.108, +0.042, 
and —0.053 mag in [/, Bw, R, z, J, //, and AT^, respec- 
tively. 

Finally, we tested the absolute photometric calibration 
of our AGES sample relative to the SDSS. We used the 
best-fitting kcorrect model to synthesize observed pho- 
tometry in the SDSS filters for each galaxy in our AGES 
sample. We then queried the SDSS/DR7 database for 
ugriz Petrosian magnitudes of the sample, and compared 
the synthesized and observed magnitudes. We found a 
negligible (< 0.01 mag) systematic difference to « ~ 18.5, 
and a median difference of 0.01 — 0.02 mag to i ~ 21, with 
a random scatter of 0.05 mag. Therefore, we conclude 
that the absolute photometric calibration of our AGES 
photometry is consistent with the SDSS to better than 

3.3. Stellar Masses 

We estimate stellar masses for the galaxies in our sam- 
ple using iSEDf it, a Bayesian SED-fitting code that uses 
population synthesis models to infer the physical prop- 
erties of a galaxy given its observed broadband SED (see 
J. Moustakas et al. 2012, in prep., for additional de- 
tails). As in §2.2.21 we adopt the BC03 population syn- 
thesis models based on the Padov a isochrones, the STELIB 
stellar library, and the iChabried ()2003D IMF. Given the 
broadband fiuxes Fi of a galaxy at redshift z in i = 1, iV 
filters, iSEDf it uses Bayes' theorem to compute the pos- 
terior probability distribution function (PDF) 



p{Q\Fi,z)=p{Q)xp{F„z\Q) 



(2) 



(iKa uffmann et al 
1200 



for a given set of model parameters Q (stellar mass, 
age, metalhcity, etc.). Here, p{Q) is the prior proba- 
bility of the model parameters, p{Fi,z\Q) is the like- 
lihood C oc exp[— x^(/i, z, Q)/2] of the data given the 
model, and is the usual goodness-of-fit statistic appro- 
priate for normally distributed photometric uncertain- 
ties. We draw each model parameter from a specified 
prior probability distribution using a Monte Carlo tech- 
nique, which is equivale nt to multiplying the likelihood 
by the prior probability (jWalcher et al.ll20lol ). Once 
has been computed for every model, the marginalized 
posterior PDF of the parameter of interest, for exam- 
ple p{M) for the stellar mass, follows from equation ([2]) 
after integrating over th e other "nuisance " parameters 
2 003al iSalim et al.ll2007l: lAuger et al.l 
The advantage of this ap- 
proach over traditional best-fitting (maximum likeli- 
hood) techniques is that it accounts for both photometric 
uncertainties and physical degeneracies among different 
models. We adopt the median of the posterior PDF as 
the best estimate of the stellar mass, and the uncertainty 
as 1/4 of the 2.3 — 97.7 percentile range of the p{M) dis- 
tribution, which would be equivalent to Icr for a Gaussian 
distribution. 

We assume exponentially declining star formation his- 
tories of the form ip{t) oc exp(— t/r), but allow each 
model galaxy to experience one or more stochastic bursts 
of star formation. We draw 1/r from a uniform dis- 
tribution in the range 0.01 — 10 Gyr~^, corresponding 
to a characteristic timescale for star formation ranging 
from r = 0.1 Gyr (similar to an instantaneous burst), 
to r = 100 Gyr (continuous star formation). We al- 



iTavlor et all I20T1I) 
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low the age t (time for the onset of star formation) 
of each model to range with equal probability between 
0.1 — 13 Gyr, although we disallow ages older than the 
age of the Universe at the redshift of each galaxy. We 
assume a uniform prior on stellar metallicity in the range 
0.004 < Z < 0.0 4, and adopt the time -dependent atten- 
uation curve of iCharlot fc Fall (|2000f ). in which stellar 
populations older than 10 Myr are attenuated by a factor 
fi times less than younger stellar populations. We draw 
II from an order four Gamma distribution that ranges 
from zero to unity centered on a typical value {fi) = 0.3 
(|Charlot fc Longhettil [200TI: IWild et al.llMTh . and the 
V-band optical depth from an order two Gamma distri- 
bution that peaks around Ay ~ 1.2 mag, with a tail 
to Ay ~ 6 mag. Finally, we allow each model galaxy 
to experience one or more random bursts of star forma- 
tion at time tb > t in each 2 Gyr time interval over the 
lifetime of the galaxy with a probability of 50%. We 
characterize each burst by a duration Ati,, drawn from a 
logarithmic distribution in the range 30 — 300 Myr, and a 
fractional stellar mass formed Ff,, whi ch we draw from a 
logarithmic distribution from 0.01 — 4 (jSalim et al.ll2007l : 
IWil d et al. 2009;). For simplicity, we assume that the 
stellar population that forms in the burst has the same 
stellar metallicity and is attenuated by the same amount 
of dust as the pre-burst stellar population. 

The largest systematic uncertainty affecting our stel- 
lar mass estimates is likely due to our choice of IMF, al- 
though if the IMF is spatially and temporarily invariant 
(i.e., universal) as is commonly assumed, then adopting 
a different IMF should simply shift all our ste l lar m asses 
by a fixed amo unt. For example, the lSalpeteii (|1955[ ) and 
iKroupal (120011) IMFs are 0.25 dex and 0.03 dex heavier 
than the IChabrie'^ (pOOl IMF, re spectively, wh i le the 
so-called "diet Salpeter" IMF of iBcll fc de Jond (f200l 
is 0.1 dex lighter. Although there are recent tantalizing 
suggesti ons of a Salpeter (or heavier) IMF in early-typ e 
galaxies (Ivan Dokkum fc Conrovll2010l:lTreu et al.lf2010l) . 
current observations heavily favor a universal Chabrier- 
or Kroupa-likc IMF in nearby star-forming disk g alax- 
ies, including the Milky Way (jBastian et al.l [20101 and 
references therein). 

Varying our prior distributions — for example, adopt- 
ing different metallicity or attenuation priors, or assum- 
ing smooth, burst-free star-formation histories — changes 
our stellar masses by < 0.1 dex in the mean, and by 
less than a factor of two for individual objects. We also 
test the effect of a,dopting different population synthesis 
models (lMarastonll200l iConrov fc Gui^lMol) and find 
comparable systematic differences in the derived stel- 
lar masses for both our SDSS and AGES samples, i.e. 
typically < 0.1 dex. We refe r the interested rea.der t o 
iKamiappan fc GawiseJ (I2007D . iMarchesini etall (|2009f ). 
and uzzin et al.l (|200^ for recent detailed discus- 
sions of the uncertainties associated with deriving stellar 
masses from broadband photometry. 

4. SAMPLE SELECTION 

The AGES main galaxy sample contains 12, 473 galax- 
ies brighter than Iab = 20.45 with we ll-measured red- 
shifts over 7.9 deg^ of the Bootes field (jKochanek et al.l 
120 lit) . We exclude from this sample five Hectospec con- 
figurations that could not be flux-calibrated (see ^2.2.1|) . 
and two additional configurations with very poor flux- 



TABLE 1 
AGES AND SDSS Galaxy Samples 



Sample 




Section'' 


AGES 








Parent 


10838 




El 


Emission-Line 


4033 




4.1 




SF'^ 


3205 




4.2 




AGN 


828 




O 




Abundance"^ 








KK04 


2975 




5.3 




T04 


3191 




5.3 




M91 


2969 









SDSS 



Parent 


427504 


2.3 


Emission-Line 


108040 


4.1 


SF'^ 


76411 


4.2 


AGN 


31629 


4.2 


Abundance'^ 






KK04 


75378 


5.3 


T04 


75976 


5.3 


M91 


75313 


5.3 



^ Number of galaxies in each sample. 
^ Section containing details regarding 
how we select each sample. 
Star-forming galaxies. 
Number of galaxies with well- 
measured oxygen a bu ndanc es based 
on the nVIcGaugbl 119911 . M91), 
ITrcmonti et al. ! Il200l. T04 ), and 
n<:obulnickv & Kewlevi l|200l . KK04) 
abundance calibrations. 

ing. However, we preserve the statistical completeness 
of the sample by upweighting the remaining objects on 
a field-by-field basis to account for the missing config- 
urations. Finally, we restrict our analysis to galaxies 
with 15.45 < Iab < 20.4 and 0.05 < z < 0.75, leav- 
ing 10,838 galaxies (Table |T]), or an effective sample of 
24, 525 galaxies after correcting for sparse-sampling and 
other sources of incompleteness (see ii2.2.1[) . In H4.ll we 
select a subset of these objects with the requisite nebular 
emission lines to estimate their oxygen abundances, and 
in M.2I we apply a variety of multiwavelength criteria to 
identify and remove AGN from our sample. Finally, in 
§4.31 we present some of the basic properties of our final 
sample of star-forming emission-line galaxies, and define 
volume-limited subsamples. 

4.1. Selecting Emission-Line Galaxies 

We select a sample of emission-line galaxies based on 
the strength of H/3. The principal advantages of using 
11/3 are that its strength is proportional to the instanta- 
neous SFR, and it is less sensitive to variations in dust 
attenuation, excitation , and metallicity than other strong 
optical emission lines f Kennicuttlll992t iMoustakas et al.l 
12006; Gilbank et al.„2010i r From a practical perspective, 
II/3 also has the advantage that it is observable in both 
AGES and the SDSS over the full redshift range of inter- 
est. 

In the top panel of Figure |2| we plot the 11/3 luminos- 
ity, L(II/3), versus redshift for all the galaxies in AGES 
with a significant H/3 emission line. For simplicity we 
ignore both relative and absolute aperture effects in our 
selection (see also §7.4p . We select all objects {small black 
points) above the black curve, which corresponds to a flux 
hmit of F[ll/3) > 3 x 10"^'' erg s~^ cm~^. In the bottom 
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Redshift 

Fig. 2. — Selection of emission-line galaxies in AGES based on the 
H/3 line-strength. In the top panel we plot H/3 luminosity, L(H/3), 
vs. redshift, and seleet all galaxies [small black points) with an 
observed H/3 line-flux -F(H/3) > 3x10"^^ erg s~^ cm"^ (solid black 
curve). The bottom panel plots rest-frame EW(H/3) vs. redshift for 
the resulting emission-line sample, where for comparison we also 
show the objects that fail our fiducial flux cut (light blue points). 
Our selection on H/9 luminosity is physically motivated because 
L(H/3) is proportional to the instantaneous SFR. 

panel we plot rest- frame EW(H/3) against redshift for the 
resulting emission-line sample, where for comparison we 
also show the objects that fail our fiducial flux cut {light 
blue points). The EW(H/3)-redshift plot shows that our 
sample spans a wide range of EW(H/?) at all redshifts, 
although the distribution narrows at higher redshift due 
to the increasing limiting L(H/3) luminosity. 

In order to estimate the gas-phase oxygen abundances 
of the galaxies in our sample we require a measurement 
of the [O ii] and [O ill] A5007 emission lines (sec g5.2p . 
In Figure Owe plot [O Iij/H^ and [O ill] /H/3 versus red- 
shift for our II/3-selected AGES sample. The small black 
points correspond to well-measured line-ratios, while the 
orange arrows represent upper limits. Focusing on [O ii] 
first (Fig.[3l top), we find that the data exhibit a distinct 
lower envelope around log([0 ii]/II/3) = —0.3 (dashed 
line), below which we find only ~ 6% of our sample, in- 
cluding nearly all the upper limits. Visual inspection of 
the spectra with low [O ii]/II/3 ratios reveals that many 
exhibit a highly reddened stellar continuum, i.e., they are 
dusty. However, because these objects constitute such a 
small fraction of the parent sample, and because their 
redshift distribution is consistent with the redshift dis- 
tribution of the parent H/3 sample, we do not expect 
that excluding them will bias our results. Therefore, 
in addition to our H/3 flux cut, we further restrict our 
sample to have a well-measured [O ii] emission line with 
log([0 II] /H/3) > -0.3. 

In contrast to [O ii]/H/3, the [O ill] /H/3 ratios for ap- 
proximately one-quarter of the H/3-selected AGES galax- 
ies are upper limits (Fig. [31 bottom). This result is not 
terribly surprising because [O III] is especially sensitive to 
variations in metallicity and excitation; in particular, the 
[O iii]/H/? ratio decreases precipitously with increasing 



meta llicity and decreasing excitation (iKewlev fc Dopital 
I2002f) . However, we show in ^7.11 that the potential 
bias against metal-rich galaxies in an [O ill] -limited sam- 
ple does not significantly affect our measurement of the 
A4 — Z relation. Therefore, to significantly simplify the 
subsequent analysis, we proceed by requiring the galax- 
ies in our emission-line sample to have a well-measured 
[O III] emission line (the black points in Fig.[3l bottom). 

The final number of emission-line galaxies in AGES 
satisfying our H/3, [O ii], and [O ill] selection criteria is 
4033, or 37% of our parent sample. For our SDSS sam- 
ple we apply the same selection criteria described above 
with the exception of a correspondingly brighter H/3 flux 
limit, F(H/3) > 1 x 10"^*^ erg s'^ cm-^, leaving 108,040 
galaxies, or 25% of the parent sample. For reference, 
not requiring our SDSS sample to have a well-measured 
[O III] line would increase the sample by approximately 
one-third. In Table [T] we summarize the total number of 
objects in these and all subsequent subsamples. 

4.2. Identifying and Removing AGN 

The Balmer and coUisionally excited forbidden lines 
present in the integrated optical spectra of galaxies can 
arise from photoionization by young, massive stars, or 
by the non-thermal photoionizing continuum of an AGN. 
Therefore, when deriving abundances from optical emis- 
sion lines it is crucial that AGN be identified and re- 
moved from the sample. We leverage the tremendous 
multiwavelength coverage of the Bootes field to identify 
AGN using a variety of criteria based on optical. X-ray, 
mid-IR. and radio observations. 

Beginning with our parent sample of emission-line 
galaxies (see §4.ip . we identify type 1 (broad- line) AGN 
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Fig. 3.— Metallicity-sensitive [O ll]/H/3 and [O III] A5007/H/3 
line-ratios vs. redshift for our H/3-limited sample of galaxies (see 
Fig-0- We represent galaxies with well-measured line-ratios using 
small black points and upper limits using orange arrows. In the 
upper panel the dashed line at log([0 ll]/H/3) = —0.3 represents 
the criterion we use to select galaxies with significant [O ll] emis- 
sion, which eliminates just ~ 6% of the H/3-selected sample. By 
contrast, the bottom panel shows that the fraction of galaxies with 
only upper limits on [O III] A5007 is not negligible due the greater 
sensitivity of [O III] to variations in metallicity and excitation. 



Evolution of the A4 — Z Relation 
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-1.5 -1.0 -0.5 0.0 0.5 
log ([N II] A6584/Ha) 

Fig. 4. — AGES emission-lino diagnostic (BPT) diagram used 
to classify our sample into star-forming galaxies or AGN based 
on their position below or above the solid curve, respectively 
UKauffmann et al.l2003g '). The contours, enclosing 25%, 50%, 75%, 
and 95% of the sample, represent galaxies with well-measured 
[N Il]/Ha and [O lll]/H/3 ratios; the small black points are individ- 
ual galaxies outside the 95% contour level. The arrows represent 
objects with la upper limits on either: (1) [N Ilj/Ho? and [O lll]/H/3 
{orange left- and down-pointing arrows); (2) just [N Il]/Ha {blue 
left-pointing arrows); or (3) just [O lll]/H/3 {green down-pointing 
arrows). Objects with upper limits on any of these line-ratios can 
be classified as star-forming if they lie in the star-forming region 
of the BPT diagram; otherwise, their classification is ambiguous. 
However, objects with upper limits on [O lll]/H/3 can be classified 
as AGN nearly unambiguously if they have log([N Il]/Ha) > —0.3 
{vertical dotted line). 

as objects with broad (typically, a > 800 km s~^ 
FWHM) Ha, H/3, or Mg ii A2800 emission, depend- 
ing on which lines are observable given the AGES 
spectral range. In our SDSS sample broad-line AGN 
have been removed already. In addition, we re- 
move any object with significant [Ne v] A3426 emis- 
sion, which is a hallmark of AGN activity, or an 
enhanced [Ne ill] A3869/II/3 ratio, which is also 
indic ative of a non-thermal photoionizing radiation 
field (|Rola et al.l [1997; Dc ssauges-Zavadsky ct al. 200(i 
IPerez-Montero et al.ll2007[) FI 

Next, we use the standard [O ill] A5007/II/3 versus 
[N ii] A6584/IIa emission-line diagnostic d iagram (Fig.Hl) 
to identify type 2 (narrow- l ine) AGN (iBaldwin et alj 
[1981 IVeiheux fc Osterbrockl 1198% IKewlev et al.l I2001t 
hereafter, the BPT diagram). We do not show the SDSS 
BPT diag ram because it exhibits t he same basic features 
(see, e.g.. iKauffmann et al.ll20033 ). We identify objects 
located below and to the left of the empirical d e marca - 
tion {solid curve) proposed bv IKauffmann et all ()2003cD 
as sta r-forming galaxi e s, and the remaining objects as 
AGN. iStasinska et al.l (|2006( ) estimate that the maxi- 
mum contribution to the integrated II/3 line-flux from 
an AGN for an obje c t classi fied as star-forming using 
the IKauffmann et all (|2003c|) criterion is ~ 3%. For 
comparison, the "maximum starburst" curve proposed 

Strong [Ne ni] A3869 emission is also observed in metal-poor 
dwarf galaxies (e.s.. Hzotov fc Thuan|[l998l) : however, these types 
of objects are only present in our sample at low redshift and easily 
identified through visual inspection of their optical spectrum. 



by IKewlev et all (I2001D allows up to a 20% - 30% con- 
tribution to H/? f rom an AGN. We adopt the stricter 
IKauffmann et al] ()2003cD demarcation to minimize the 
metallicity bias incurred from contamination to the neb- 
ular emission lines from AGN (see ^7.3p . 

In Figure Ul we differentiate between four different sce- 
narios: (1) galaxies with well-measured [N ii]/IIa and 
[O iii]/II/3 line-ratios {contours and greyscale); (2) ob- 
jects with [N ii]/Ha upper limits {blue arrows); (3) ob- 
jects with [O lll]/H/3 upper limits {green arrows); and (4) 
galaxies with both [N Iij/Ha and [O iii]/H/3 upper limits 
{orange arrows). An object with an upper limit on any of 
these line-ratios can be classified as star-forming if it lies 
in the star-forming region of the BPT diagram; other- 
wise, its classification is ambiguous, with one exception: 
a galaxy can be classified as an AGN if it has an upper 
limit on [O iii]/H/3 but a well- measured [N ii]/IIa ratio 
with log([N ii]/Ha) > —0.3 {vertical dotted line). 

In order to classify the rest of the AGES sample, in- 
cluding all the objects at z > 0.4 (which lack [N ii]/IIa 
measurements), we must rely on other diagnostics of 
AGN activity. Recognizing the difficulty of measur- 
ing the [N ii]/Ha ra tio of intermediate-redshift galax- 
ies, lYan et all (|2011[ ) have proposed a new BPT-like 
diagnostic diagram that replaces the [N ii] /Ha ratio 
with the rest-frame U — B color (see iLamareilld [2010l 
iLara-Lopez et al.ll2010al and lJuneau et al.ll2011l for alter- 
native empirical diagnostic diagrams) . The physical mo- 
tivation behind this classification scheme is that AGN are 
more likely to be found i n galaxies dominated by evolved 
stellar populations (e.g.. IKauffmann et al.|[2003c|) : there- 
fore, at fixed [O iii]/H/3 ratio, galaxies whose integrated 
emission-line spectrum is dominated by star-formation 
are more likely to be blue in their U — B color. 

In Figure [5] we plot U — B versus [O iii]/H/3 for our 
{left) SDSS and {right) AGES galaxy samples. The light 
blue contours represent objects classified as star-forming 
galaxies using the BPT diagram, the grey contours repre- 
sent objects classified as AGN, and the dark blue points 
indicate objects that cannot be classified using the BPT 
diagram. For clarity we do not include in this figure the 
small number of SDSS galaxies that could not be classi- 
fied using the BPT diagram, although nearly all of them 
would be classified as AGN using the lYan et afl ((20lTI) 
diagnostic diagram. Among our AGES galaxies, the tan 
squares represent objects identified as AGN using vari- 
ous other optical. X-ray, mid-IR, and radio criteria, as 
described below. Fin ally, the so l id line is the empirical 
boundary proposed bv lYan et al.l ()2011t) to separate star- 
forming galaxies from AGN, which are located below and 
above the solid line, respectively. 

Examining Figure [3 we find that the lYan et all (|2011l ) 
diagnostic diagram successfully classifies as star-forming 
more than 98% of the BPT-classified star-forming galax- 
ies; however, it also identifies as star-forming roughly 
50% of the objects classified as AGN using the BPT 
diagram. In other words, the diagram has high com- 
pleteness for star-forming galaxies, but at the expense 
of a relatively high contami nation rate (i.e., l ow purity). 
However, as emphasized bv lYan et al.l (|2011l ). nearly all 
of these "misclassified" objects are, in fact, composite 
systems with an admixture of star-formation and AGN 
activity. Moreover, the contribution from the AGN to 
the integrated emission-line spectrum in composite sys- 
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Fig. 5. — BPT -like diagnostic diagram for our (left) SDSS and (right) AGES galaxy samples. This diagram has been proposed by 
lYan et ahl | |2011I '| as an empirical means of segregating star- forming galaxies {below the solid line) from AGN {above the solid line) when the 
ionization-sensitive [N Il]/Ha ratio cannot be measu red. The blue contours enclose 25%, 50%, 75%, and 95% of the galaxies classified as 
star-forming using the BPT diagram (see Fig.|4]and i|4.2ll . while the grey contours show the distribution of BPT-classified AGN. The dark 
blue points represent galaxies that cannot be classified using the BPT diagram. For clarity we do not plot the small number of unclassified 
SDSS galaxies, although nearly all of them lie in the AGN portion of this diagnostic diagram. Within the AGES sample, the tan squares 
represent galaxies classified as AGN using other multiwavclength (X-ray, mid-IR, and radio) criteria. With respect to objects classified as 
star-forming galaxies using the BPT diagram, we find that this technique has high completeness, but low purity. In other words, while 
> 98% of BPT-classified star-forming galaxies are su ccessfully c l assifie d as star-forming using this diagnostic diagram, more than 50% of 
AGN are also identified as star-forming. As noted bv lYan et al.l l)20111 V however, the bulk of these objects are composite systems, with an 
admixture of ongong star-formation and AGN activity. 



terns is always sub-dominant (5% — 30%: iStasinska et aTl 
[2006h . 

To mitigate the tendency for the lYan et all ()2011f ) 
diagnostic diagram to classify composite systems as 
pure star-forming galaxies, we utilize a variety of com- 
plementary multiwavclength criteri a to identify AGN 
in our AGES sample (see also Hickox et al.l 120091 : 
iRuiopakarn et al.1 l2010t lAssef et all 1201 ID . First, we 
classify as an AGN any object with more than a 50% 
probability of being an X -ray point source (jKenter et al.l 
120051 iBrand et al]l2006[ ). Although the X-ray emission 
in some of these objects may aris e from processes as - 
sociated with star formation (e.g.. iWatson et aril2009[ ). 
roughly 90% of them have total 0.5 — 7 keV X-ray lumi- 



nosities Lx ^, 



10^^_erg_s ^ 



consistent with hosting an 



AGN (iGcorgakakis et al.l 120061 ). We also isolate objects 
lying in the "Stern wedge" of the IR AC [3.6] - f4 - 5l ver - 
sus [5.8] - [8.0] color-color diagram (jStern et al.ll2005f ). 



scured and unobscured in the X-ray 1 


Hickox et al.1120071: 


iGoriian et al.ll2008HAshbv et al.ll200S 


) . Finally, we com- 



pute the radio power at 1.4 GHz, Pi. 4, of the radio 
sources in our sample (see §2.1.2p . if-corrected assuming 
an intrinsic power-law spectrum Si, oc with a = 0.5 
(jPrandoni et al.l[2006l ). and iden tify AGN as sources hav- 
ing log P^ A > 23 .8 W Hz"! (|Kauffmann et al.l [200I 
IHickox et al1[2009l) . 

One potential objection to the preceding analysis is 
that AGN identified at X-ray, mid-IR, and radio wave- 
lengths should not necessarily be removed from an optical 
emission-line abundance study. For example, there exists 



an interesting (albeit rare) class of objects with signifi- 
cant X-ray emission but with no discernable evidence of 
AGN activity in the optical, so-call ed X-ray Bri ght Opti- 
cally Normal Galaxies (XBQNGs: iComastri et al. 200i 
iCivano et all [20071 : lYan et al.l l20m Similarly, AGN 
identified based on their mid-IR colors may be so ob- 
scured at optical wavelengths as to not affect the opti- 
cal emission lines in any important way. Nevertheless, 
we conservatively reject any object (not otherwise classi- 
fied as star-forming using the BPT diagram) whose op- 
tical emission lines might be contaminated by an AGN. 
The reason for adopting this position is that the oxygen 
emission-line ratios of unidentified AGN are similar to 
those of metal-poor star-forming galaxies, and therefore 
could bias the amount of metallicity evolution that we 
infer, a point that we revisit in i )7.3l 

To summarize, our final SDSS and AGES star-forming 
galaxy samples contain 76, 411 and 3205 galaxies, respec- 
tively (see Table H]). However, 432 (14%) of the AGES 
emission-line galaxies, 85% of whic h are at z > 0.4, w ere 
classified as star-forming using the lYan et al.l (|2011[ ) di- 
agnostic diagram, which carries a potential ~ 50% con- 
tamination rate from composite (i.e., low- level) AGN. 
On the other hand, by using additional multiwavelength 
AGN diagnostics it is likely that the contamination rate 
in our sample is much lower. Nevertheless, we conserva- 
tively estimate that < 10% of the galaxies in our star- 
forming galaxy sample, or up to ~ 50% of the objects 
at z > 0.4, may host weak (composite) AGN. We show 
in §7.31 however, that these unidentified optical AGN, if 
present, have a minimal effect on our conclusions. 
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4.3. Sample Properties 

In Figure H] we compare the distributions of redshift, 
absolute _B-band magnitude, U — B eolor, and stellar 
mass for the star-forming galaxies in AGES {filled light 
blue histograms) against the full parent sample [unfilled 
light grey histograms)^ both corrected for the statisti- 
cal incompleteness of the survey (see i 32.2.1|) . We also 
show for comparison the distribution of properties for 
our SDSS star-forming galaxy sample (unfilled dark red 
histograms), where we have scaled the SDSS distribu- 
tions down by a factor of 15 for visualization purposes. 
Not unexpectedly, our emission-line cuts preferentially 
remove massive, luminous galaxies with red optical col- 
ors, the majority of which are presumably old, early-type 
galaxies with little or no ongoing star formation (e.g., 
iKauffmann et al.ll2003bl : lZhu et al.ll2010L but see the dis- 
cussion in ^7.ip . Corrected for incompleteness, the me- 
dian (mean) redshift of our SDSS sample is Zmed = 0.087 
(0.094), compared to z^cd = 0.27 (0.30) for our AGES 
sample; otherwise, the distributions of luminosity, stellar 
mass, and color for the two samples are very similar. 

In Figure [7] we examine the properties of the star- 
forming galaxies in AGES in more detail by plotting B- 
band absolute magnitude and stellar mass versus red- 
shift. The distribution of points follows the charac- 
teristic shape of a flux-limited survey, in which intrin- 
sically luminous, massive galaxies are observable over 
the full range of redshifts, unlike low-luminosity galaxies 
(i.e., Malmquist bias). This figure also reveals the rel- 
ative dearth of luminous, massive star-forming galaxies 
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Fig. 6. — Distribution of redshift, absolute B-band magnitude, 
rest-frame U — B color, and stellar mass for the AGES parent sam- 
ple {unfilled light grey histograms), and our sample of star-forming 
emission-line galaxies [filled light blue histograms), weighted to cor - 
rect for the statistical incompleteness of the survey (see ^2.2.11 1. 
The unfilled dark red histograms show for comparison the distri- 
butions of the same properties for our SDSS star-forming galaxy 
sample, where we have scaled the SDSS distributions down by a fac- 
tor of 15 for comparison purposes. As expected, our emission-line 
selection criteria preferentially remove massive, luminous, optically 
red galaxies. 
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Fig. 7. — Absolute B-band magnitude and stellar mass vs. red- 
shift for the star- forming galaxies in AGES. The black boundaries 
demarcate volume-limited subsamplcs in si x red shift bins centered 
on 2 = 0.1, 0.2, 0.3, 0.4, 0.5, and 0.65 (see [gjland Tablc[2ll. 

in AGES at low redshift due to the limited cosmological 
volume probed. 

We use Figure [7] to calculate the minimum absolute 
magnitude and stellar mass above which our sample is 
volume-limited function of redshift. The limiting 
Mb is directly proportional to the survey flux limit, 
modulo small variations in iT-corrections among differ- 
ent galaxy types. However, the limiting stellar mass 
has an additional dependence on the mass-to-light ra- 
tio (i.e., star- formation history) of each galaxy. At fixed 
stellar mass and redshift a galaxy with a large (small) 
mass-to-light ratio will lie preferentially below (above) 
the survey flux limit. Since observations suggest that 
galaxies of a given m ass with high star formation rates 
tend to be metal-poor (jMannucci et al.|[2010l : lYates et"al] 
1201 in . properly accounting for the stellar mass limit en- 
sures that we do not bias our measurement of the Ai — Z 
relation toward lower metallicity. 

To calculate the absolute magnitude and stellar 
mass limits we adopt an empirical approach (see, e.g., 
ILamareille et all 120091 ). First, we divide AGES into six 
redshift bins centered on z = 0.1, 0.2, 0.3, 0.4, 0.5, and 
0.65; corresponding to approximately equal intervals of 
cosmic time (^ 1 Gyr). Next, we compute the absolute 
magnitude and stellar mass each galaxy would have if its 
apparent magnitude, /, equalled the magnitude limit of 
the survey, /AB.iim = 20.4. Specifically, we calculate 

A^S,mi„ = Mb -{I- /lim) (3) 

and 

log (Mmin/MQ) = log (M/Mq) + 0.4 (/ - hm) (4) 

for each galaxy. The distributions of MB,min and 
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TABLE 2 

AGES Limiting _B-Band Magnitude & 
Stellar Mass'' 



Redshift Range 






0.05 - 0.15 


-17.47 


8.75 


0.15 - 0.25 


-18.60 


9.38 


0.25 - 0.35 


-19.58 


9.88 


0.35 - 0.45 


-20.40 


10.26 


0.45 - 0.55 


-21.06 


10.52 


0.55 - 0.75 


-21.76 


10.67 



* Absolute B-band magnitude and stellar mass 
as a function of redshift above which the 
AGES star-forming galaxy sample is more than 
50% complete, accounting for variations in K- 
corrections and mass-to-light ratio. 

■Mmin in each redshift shcc reflect the variations in 
/^-corrections and mass-to-hght ratio among different 
galaxies. We then define MB,iim and Mum as the lu- 
minosity and stellar mass that includes more than 50% 
of the galaxies. In Figure [7] we show the boundaries that 
define our volume-limited samples in each redshift inter- 
val, and list the 50% completeness limits in Tabled] For 
reference, the 75% and 95% stellar mass completeness 
limits are roughly 0.15 dex and 0.4 dex higher, while 
the corresponding limiting Mb values are approximately 
0.1 mag and 0.35 mag brighter, respectively. 

5. DERIVING GAS-PHASE OXYGEN 
ABUNDANCES 

5.1. Background on Strong-Line Metallicity 
Calibrations 

The observational technique of using optical emission 
lines to derive the physical conditions (density, temper- 
ature, chemical composition) of the ionized gas in H ii 
regio ns and star - forming gala,xies has a long, rich his- 
tory (|Alledll942t IShielddllQM IStasinskallMm . In par- 
ticular, considerable effort has gone into using optical 
spectroscopy to infer the abundance of oxygen in star- 
forming regions. Oxygen is important because it is the 
most abundant clement (b y number) in the Un iverse af- 
ter hydrogen and helium (jAsplund et al.|[2009f ). making 
it a crucial coolant in the interstellar medium; moreover, 
all the principal oxygen transitions are observable in the 
rest-frame optical. 

The classical, so-called direct method of inferring the 
gas-phase oxygen abundance is to measure the elec- 
tron temperature of the ionized gas from the strength 
of an auroral line, usually [O ill] A 4363, relat i ve to 
a nebular line such as [O ill] A5007 ("Dincr steinl 119901: 
IGarnettl[200l lOsterbrock fc Ferland 2006)0 Unfortu- 
nately. [O III] A4363 and other temperature-sensitive au- 
roral lines (e.g., [N ii] A5755 [S ill] A6312, [O ii] A7325) 
become vanishingly weak with increasing metallicity; 
therefore, they are rarely detected in the integrated spec- 
tra of galaxies observed as part of flux-limited surveys 
such as AGES and the SDSS, which typically include rel- 

■^^ The abundance of oxygen and other heavy elements can be 
inferred from the optical metal-recombination lines as well. How- 
ever, because these lines are extraordinarily weak, they have been 
measured in only a relatively small number of Gala ctic and extra- 
galac tic H II regions and starburst dwarf galaxies IIEsteban et al.l 
12003 . and references therein). 



atively luminous, metal-rich galaxies. For example, elec- 
tron temperatures have been measured for fewer than ~ 
0.1% of galaxies in the entire SDSS databa se, all of which 
are relatively nearby, m etal-poor dwarfs (|Kniazev et al.l 
[200l llzotov eral] [200l. 

In light of the observational challenges confronting 
direct, electron-temperature abundance measurements, 
complementary methods have been developed that rely 
exclusively on the strongest emission lines in the rest- 
frame optical. These so-called strong-line methods 
are statistical correlations between various line-ratios 
and the gas-phase oxygen abundance, 12 -|- log(0/H). 
Unfortunately, different calibrations yield widely dif- 
ferent oxygen abundances based on the same input 
emission-line spectrum. The magnitude of the system- 
atic differences in metallicity depends on which cali- 
brations are being compared, but they can be as large 
as ~ 0.7 dex, or a factor of ^ 5 (iKennicutt et all 
200l [Pirei-lVIontero fc Dia3 [20&1 iLiang et al I l2006bl: 



Nagao et all 120061: iB rcsolin et al.l [2004 120051: iBresolin 
2007t iKewlev fc Ellis on 200 ^ iLopez-Sanchez fc Esteban 
2OI0I : iStasihskal I2010t iMoustakas et al.ll2010[ ). In par- 
ticular, so-called theoretical strong-line methods, which 
are based on ab initio photoionization model calcula- 
tions fe.g.. lMcGaugh|[l99lt iCharlot fc Longhettil [2OOII: 
IKewlev fc Dopital |2002[ ) . yield oxygen abundances that 
are, on average, several times higher than abundances 
derived using either the direct method, or empirical 
strong-line methods, which are calibrated aga inst H 11 
regions with direct abundance estimates (e.g., iPilvuginI 
2OOOI I2OOII: iDenicolo et all I2OO2I: iPettini fc Pagell 1 20041: 



Perez-Monter o fc Dfad 120051: 
Stasiiiska 2006: lYin et all 120071 ) 



Pilvugin fc ThuanI 200 



5.2. Deriving Metallicities Using the R23 Parameter 

Because of the poorly understood systematic uncer- 
tainties in the nebular abundance scale described in 
it is important that the requisite emission lines we use to 
derive metallicities are observable over the full redshift 
range of interest, z = 0.05 — 0.75. The combination of 
nebular lines that satisfies this requirement f or our SDSS 
and AGES samples is the R23 parameter (|Pagel et al.l 

IMi), 



R23 — 



[O II] A3727 + [O III] AA4959, 5007 
up A4861 ■ 



(5) 



The i?23 parameter is a popular and widely used abun- 
dance diagnostic for several reasons. First, its strength is 
directly proportional to both principal ionization states 
of oxygen, and therefore does not require large, uncer- 
tain ionization corrections; the relative amount of neu- 
tral oxygen and O"^"*" is negligible in star-forming re- 
gions across a wide range of physical conditions. Second, 
i?23 is explicitly related to the oxygen abundance, un- 
like other popular diagnostics that depend implicitly on 
the relative abundance of additional elements like nitro- 
gen (see ^7.2[) . Finally, R23 can be measured from the 
ground using a combination of optical and near-infrared 
spectroscopy from z = — 4, which allows the chemi- 
cal enrichment histories of galaxies to be studied using 
a consistent met allicity diagnostic over 90% of t he age 
of the Univ erse ('K cnnicuttlll998al: iMaiolino et"all 120081: 
iRichard et al...2011. ). 
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Fig. 8. — Emission-line diagnostic di agra m used to assign galaxies to the upper, log ([N Il]/Ha) > —1.1, and lower log ([N Il]/Ha) < —1.1, 
branches of the iJ23-0/H relation (see ^5.21 for our (left) SDSS and {right) AGES star-forming galaxy samples with wcll-mcasured [N iiJ/Hq 
ratios. For reference, the contours enclose 25%, 50%, 90%, and 97.5% of the sample, and the small grey points represent galaxies that 
lie outside the 97.5% contour level. This diagram demonstrates that the galaxies in both samples overwhelming belong to the upper i?23 
branch. 



To explore the effect on our conclusions of choosing a 
particular strong-line calibration, we carry out our anal- 
ysis using three independent the o retical strong- line cali- 
bratio ns of feci: iMcGaughl (ll99lD [Kobulnickv fc KewlevI 
(j2004D . and iTremonti et all (j2004[) . hereafter referred to 
as the M91, KK04, and T04 abundance calibrations, re- 
spectively. At a fixed value of R23, both observations and 
photoionization modeling show that there is a second- 
order dependence of 0/H on the ionization par ameter 



U, or the hardness of the ionizing radiation field (IShields 
199n[ IMcGaughl [T991I : IPilvuginI 120011 : iKewlev fc Dot^ 



20021) . The M91 and KK04 calibrations constrain the 
ionization parameter using the O32 ratio, where 



032 = 



_ [O III] AA4959, 5007 



[O II] A3727 



(6) 



while the T04 calibration ignores this effect. All three 
calibrations are based on photoionization models, and 
therefore they yield metallicities that are considerably 
higher than those implied by the d irect method or empir- 
ical strong-line cal i bratio ns (e. g.. iPettini fc Page| 120041 : 
IPilvugin fc Tbiaa [2005t see IKewlev fc EUisonl \20m . 
However, as iMoustakas et al. 1 mM 

argue, existing em- 
pirical strong-line calibrations should be applied to in- 
tegrated galaxy spectra with great care because they 
are being extrapolated into a part of physical parame- 
ter space that is not well-constrained by current obser- 
vations of H II regions. In Appendix El we compare the 
M91, KK04, and T04 calibrations in more detail. 

Despite its popularity and wide-spread use, i?23-based 
abundances suffer from several potential systematic un- 
certainties. First, the relationship between R23 and 0/H 
is famously double- valued (see, e.g., Fig. IA1[) : a given 
value of i?23 can correspond to a solution on the metal- 
rich upper branch, or on the metal-poor lower branch; the 
region where the two branches overlap is known as the 



turn-around region. Fortunately, the bulk of the galaxies 
in our SDSS and AGES samples are relatively luminous 
and massive, and therefore overwhelmingly belong on the 
upper branch. We illustrate this point in Figure [HI where 
we plot [N II] A6584/Ha versus i?23 for our (left) SDSS 
and (right) AGES star-forming galaxy samples. Ob- 
jects that belong on the lower i?23 branch generally have 
log ([N ii]/Hq;) < —1.1, w hile upper-branch galaxies have 
log([N ii]/Hq) > -1.1 (IMcGaughl [I991I IContini et al.l 
[200a IKewlev fc EUisonl l2008| ). We find that fewer than 
1% of SDSS galaxies, and only ~ 3% of AGES galaxies 
belong to the lower R23 branch. Of course, we do not 
have a measurement of [N 11] /Ha for galaxies in AGES 
at 2 > 0.4 because of the red wavelength cut-off of our 
spectra. However, we assume that all these galaxies lie 
on the upper i?23 branch, which is not an unreasonable 
assumption given their luminosities and stellar masses 
(sec Fig. [7]). In particular, this assumption is conserva- 
tive because it minimizes the amount of chemical evolu- 
tion these objects would have to undergo to lie on the 
present-day Ai — Z relation (see ^6.1.ip . 

The second issue to consider is that both i?23 and 
O32 are sensitive to dust attenuation. Unfortunately, we 
cannot correct for dust reddening because the Balmer 
decrement, Hq;/H/3, can be measured only for galaxies 
in AGES at z < 0.4, while the H/3/H7 ratio is gener- 
ally too noisy to give a r eliable estimate of the r edden - 
ing. Therefore, following iKobulnickv fc Phillipj ()2003D . 
we use EWs to estimate oxygen abundances. Besides be- 
ing relatively insensitive to dust extinction as we show 
below, EWs have the added advantage of being impervi- 
ous to flux-calibration issues (sec H2.2.2p . 

With some simple assumptions it is easy to show that 
the EW-bascd i?23 and O32 parameters are 



EW(i?23) 



aEW([0 II]) + EW([0 III]) 
EW(H/3) 



(7) 
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Fig. 9. — (Top) Correlation between oxygen abundances esti- 
mated using reddening-corrected flu xes a nd oxygen abundances 
derived using EWs and o = 1.0 (see ^5.2|l . This comparison uses 
the KK04 abundance calibration, but the same conclusion holds if 
we use either the M91 or T04 calibrations. (Bottom) Abundance 
residuals vs. rest-frame U — B color. The open squares give the 
median residual metallicity in 0.05 mag wide bins of U — B color 
and the error bars indicate the interquartile range. The median 
residuals vary by < 0.05 dex with a scatter of just 0.01 — 0.06 dex; 
therefore, abundances derived using EWs provide a reliable way of 
estimating the metallicity when reddening-corrected fluxes are un- 
available. For reference, the contours in both panels enclose 50%, 
75%, 90%, and 97.5% of the sample, and the solid line shows the 
one-to-one relation; the small grey points are individual galaxies 
lying outside the 97.5% contour level. 



and 



where 



aEW([0 II]) ^ 



1,(4934) 
/c(3727) 



y)[fe(4934)-fc(3727)] 



(8) 



(9) 



is a constant for each galaxy that depends on its re- 
cent star-formation history and level of dust attenua- 
tion. In equation (O, 7,(4934) and 7,(3727) are the in- 
trinsic (unattenuated) flux of the stellar continuum at 
A = 4934 A (the wavelength midway between H/3 and 
[O III] A5007) and 3727 A, respectively; AE{B-V) is the 
reddening difference between the ion ized gas and the con- 
tinmmi (in general, AE{B-V) > 0: lCalzetti et al.lll994l: 
[Calzctti 2001); and fc(4934) and fc(3727) is the value of 
the assume d attenuation law at 4934 A and 3727 A, re- 
spectively (iKobulnickv fc Phillip j [20031) . Using popula- 
tion synthesis models a nd some basic assumpt i ons ab out 
AE{B - V) and fc(A), IKobulnickv fc PhillipsI (12001 es- 
timate a = 0.84 ± 0.3. 



We obtain an independent, empirical estimate of a by 
correlating the reddening-corrected O32 ratio, (O32)cor, 
and EW(032) using our SDSS sample. First, we esti- 
mate the nebular reddening for each galaxy using the 
observed Balmer decrement and assume an intrinsic 
Hq /H/? ratio o f 2.86 . Next, we derive (032)001 using 
the 10 DonnelH (Il994f) Milky Way extinction curve an d 
Rv = 3.1 (jMoustakas et al.ll2006l:iKennicutt et al.ll2009[) . 
Equating EW(032) and (O32)cor, we find, on average, 
a = 1.0 ± 0.16. To a very good approximation, in other 
words, equations ([7]) and ([S]) with awl reduce to equa- 
tions ^ and ^ but with EWs in place of reddening- 
corrected fluxes. 

In Figure [9] we compare the EW-bascd abun- 
dances against the abundances derived using reddening- 
corrected fluxes. For this comparison we use the KK04 
calibration, although we reach identical conclusions using 
either the T04 or M91 calibration. In the bottom panel of 
Figureiniwe correlate the abundance residuals against the 
rest-frame U — B color. We find an excellent correlation 
between the two abundances, with a negligible system- 
atic difference (< 0.05 dex) and a la scatter that ranges 
from 0.01-0.03 dex at high metallicity, to 0.04-0.06 dex 
for more metal-poor galaxies, independent of the adopted 
abundance calibration. Therefore, we conclude that EWs 
(assuming a — 1) provide an effective means of estimat- 
ing the gas-phase abundances of galaxies in the absence 
of a measurement of the Balmer decrement (but for an 
alternative conclusion see ILiang et al.l[2007h . In 
we further show that using EWs in place of reddening- 
corrected fluxes has a negligible effect on the inferred 
shape and normalization of the M ~ Z and L — Z rela- 
tions. 

5.3. From Equivalent Widths to Metallicity 

Given the measured values of EW(7J23) and EW(032), 
and the inferred or assumed 7^23 branch, we derive the 
gas-phase oxygen abundances of all the star-forming 
galaxies in our SDSS and AGES samples using the 
M91, KK04, and T04 strong-line calibrations (see Ap- 
pendix |A]). To derive the uncertainty in 12 -|- log (0/H), 
we apply the Monte C arlo-based method developed by 
iMoustakas et al.l (|2010f ). which uses the statistical un- 
certainties on EW([0 11] ), EW([0 ill]) and EW(H/3) to 
obtain a realistic estimate of the abundance error for each 
object. This algorithm also relies on quantitative criteria 
to identify galaxies that lie off the adopted abundance 
calibration, typically because they have 7^23 ^ 10 (see 
Fig EH). The final number of AGES and SDSS galaxies 
with well-measured oxygen abundances is ~ 3000 and 

75, 000, respectively, where the exact number depends 
on which strong-line calibration is being used (see Ta- 
ble H]). 

6. ANALYSIS 

Armed with stellar masses, optical luminosities, and 
gas-phase oxygen abundances for both our SDSS and 
AGES samples, we now proceed to study the change in 
the chemical abundances of star-forming galaxies since 
z — 0.75. Although we focus our analysis on the evolu- 
tion of the M — Z relation, we also examine the evolution 
of the B-band L — Z relation. In fj6T] we establish the 
M. — Z and L — Z relations at z ~ 0.1 using our SDSS 
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TABLE 3 
SDSS M~ Z AND L~ Z Relations 







M- Z^ 








Calibration^ 


12 + log (O/H)* 


log(A4*/10''>!G) 7 '7'' 


CO 


Cl 


ad 


Reddening-Corrected Fluxes'' 


KKU4 


9.115 


2.043 1.41 0.07 


9.066 


-0.231 


0.16 


M91 


9.004 


3.117 1.27 0.09 


8.919 


-0.256 


0.18 


T04 


9.098 


3.475 1.22 0.10 


9.014 


-0.291 


0.20 


MPA-JHU 


9.137 


3.436 0.89 0.09 


9.026 


-0.248 


0.16 


Equivalent Widths" 


KKU4 


9.108 


2.001 1.25 0.07 


9.054 


-0.229 


0.16 


M91 


9.006 


3.151 1.06 0.09 


8.903 


-0.248 


0.17 


T04 


9.097 


3.575 1.05 0.10 


8.997 


-0.286 


0.20 



^ Results derived usi ng the IKobulnickv fc Kewle-d ll200l . KK04), IMcGauehl l lT99ll . M91), and 
ITremonti et al.l l|2004l . T04) calibratio ns o f the H23 parameter, using both reddening-corrected 
line-fluxes and equivalent widths (sec ^5.21 for details). For comparison, wc also derive the M — Z 
and L — Z relations using the oxygen abundances published by the MPA-JHU team (see i|6.1.1l l. 
^ M- Z relation given by 12 -|- log (O/H) = 12 -|- log (O/H)* - log [1 + (M* /W^ Mq)"']. 

B-band L - Z relation given by 12 -|- log (O/H) = cq + ci{Mb + 20.5). 

Residual lu scatter about the best-fitting relation in dex. 

sample, and in %.2\ we use AGES to measure their evo- 
lution from z = 0.05 — 0.75. We conclude in ^ 36. 31 by 
comparing our results with previous determinations of 
the Ai — Z relation at intermediate redshift. 



6.1. SDSS M — Z and B-Band L — Z Relations 

We begin by determining the — Z relation for star- 
forming galaxies at z '--^ 0.1 in !j6.1.11 and the L — Z 
relation in i j6.1.2l In particular, we examine the effect 
of adopting a specific abundance calibration on the re- 
sults. In both sections we derive metallicities using the 
reddening-corrected emission- line fluxes (see i 35.2|) . al- 
though in ij6.1.3l we examine the effect of using EWs on 
our results. 

6.1.1. M ~ Z Relation at z ^ Q.l 

In Figure [10] we plot the SDSS M - Z relation us- 
ing four different abundance calibrations. The first three 
panels show the M. — Z relations based on the KK04, 
M91, and T04 abundance calibrations, and the lower- 
right panel shows, for comparison, the — Z relation 
using the oxygen abundances derived by the MPA-JHU 
tea mH The MPA-JHU abunda nces are based on fitting 
the iCharlot fc Longhettil ()2001l ) photoionization models 
to all the available nebular emission lines, and therefore 
they serve as a useful cross-check of the i?23-based abun- 
da nces. We point out t hat the T04 calibration of R23, 
bv ITremonti et al.l (|2004[ ) was derived by correlating i?23 
against the MPA-JHU metallicities; therefore, the statis- 
tical agreement between the T04 and MPA-JHU M - Z 
relations is not unexpected. 

We characterize the shape of the M. — Z relation by 
computing the weighted mean metallicity in 0.05 dex 
wide bins of stellar mass, and plot the results in Fig- 
ure [To] using red circles, blue squares, green triangles, 
and black diamonds for each of the four calibrations, as 
indicated in the figure. To fit the M. — Z relation wc 
adopt the following functional form: 

^° |http : //www ■ mpa-garching ■ mpg . de/SDSS/DR7 1 




9 10 11 9 10 11 

log (M / h,t Ms) log (M / h™^ Ms) 



Fig. 10. — SDSS M — Z relation based on four different methods 
of deriving oxygen abundances from the reddening-corrected line- 
fluxes. The first three panels show the M — Z relation based on 
the three iJ23-based strong-line calibrations adopted in this paper 
(KK04, M91, and T04), while the M - Z relation in the lower- 
right panel uses the metallicit ies fo r SDSS/DR7 galaxies derived 
by the MPA-JHU team (see i|6.1.1l l. The contours enclose 25%, 
50%, 75%, 90%, and 97.5% of the sample, and the small grey 
points are individual galaxies that lie outside the 97.5% contour 
level. The symbols in each panel (as indicated in the legend) show 
the weighted mean metallicity in 0.05-dex wide bins of stellar mass, 
and the correspond ing curves show the best-fitting JV[ — Z relation 
given by equation (IIOI I. For reference, we reproduce in each panel 
the MPA-JHU M — Z relation {solid black curve). In the lower- 
right panel wc also show for comparison the origin al SDSS M — Z 
relation (based on the SDSS/DR4) published bv ITremonti et al.l 
II2004I ) {dashed purple curve). This figure demonstrates that both 
the normalization and the shape of the M — Z relation depend on 
the method used to derive oxygen abundances. 
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/ M* 

12+log(0/H) = 12+log(0/H)*-log ^1 + [y^^ 

°(10) 

where 12 + log(0/H)* and M* are the characteristic 
metalhcity and stellar mass, respectively. We discuss and 
justify this particular parameterization of the M — Z re- 
lation compared to the more commonly used polynomial 
model in Appendix [Bj Wc show the best-fitting KK04, 
M91, T04, and MPA-JHU M-Z relations in Figure [TO] 
using long-dashcd, dot-dashed, triplc-dot-dashcd, and 
solid lines, respectively, and list the best-fitting coeffi- 
cients in Table [3l Note that the formal statistical un- 
certainties on the coefficients are negligible (< 0.01 dex), 
and therefore have not been tabulated. For reference, 
we reproduce the MPA-JHU Ai — Z relation in every 
panel of Figure 1101 In addition, in the l ower-right panel 
we sh ow the M — Z relation derived by iTremonti et al.l 
()2004[) from a sample of ~ 53,000 galaxies from the 
SDSS Data Release 4 (DR4) as a short-dashed pur- 
ple line. The origin of the discrepancy between the 
M — Z relation published by [Trcm onti et al. (200^ us- 
ing DR4 and the MPA-JHU M - Z relation we derive 
using DR7 is not clear, given that the metallicities were 
deri ved using the same maxim um likelihood technique 
and lCharlot fc Longhettil ()200lD photoionization models. 
Nevertheless, evolutionary studies and t heoretical mod- 
els tha t are calibrated against the original fTremonti et al.l 
()2004t ) j\4 — Z relation should bear these differences in 
mind. 

Returning to Figure lTOl we find that all four abundance 
calibrations yield a tight, well-defined M — Z relation — 
the Icr scatter ranges from 0.07 — 0.10 dex — that in- 
creases monotonically with stellar mass and then fiattens 
by varying amounts above a characteristic mass around 
~ 10^° Mq. However, both the normalization and shape 
of the Ai — Z relation vary significantly among the dif- 
ferent calibrations, a result that has been explored in de- 
tail by Kewley & EUison (200^^. For example, the KK04 
M — Z relation is much flatter (7 is larger) , and the char- 
acteristic mass is a factor of ~ 1.6 smaller compared to 
the other three M — Z relations, although the character- 
istic metallicities of the four M — Z relations are within 
~ 0.1 dex (25%) of one another. 

Understanding the origin of the differences between the 
M — Z relations implied by these strong-line calibrations 
is beyond the scope of th is paper (but see the discussion 
in lMoustakas et aLll2010D . The main points to take away 
from the preceding discussion and Figure [TO] are that 
both the normalization and shape of the A4 ^ Z relation 
depend on the method used to derive gas-phase abun- 
dances. This result has important implications for theo- 
retical semianalytic and hydrodynamic models of galaxy 
formation that use the M — Z relation to calibrate and 
tune so-called "recipes" for star formation efficiency, su- 
pernova feedback, and g alactic winds (e.g.. iBrooks et al. 
20071: iFinlator & Davd [2008]; iDutton fc van den Bosch! 
20091 Ibavc ct al., 201 lb}) . From an observational stand- 
point, these results emphasize the importance of focusing 
on the differential metalhcity evolution, and of consider- 
ing multiple abundance diagnostics to ascertain whether 
the rate of metalhcity evolution depends on the adopted 
abundance calibration (see especially ij6.2.ip . an issue 



that heretofore has largely been ignored. 

We conclude this section by verifying that the Ai — Z 
relations measured from our SDSS sample at z ^ 0.1 and 
the subset of AGES galaxies in our lowest redshift bin, 
0.05 < z < 0.20 (hereafter, the low-redshift AGES sam- 
ple), are consistent. In FigurefTDwe plot the KK04-bascd 
Ai — Z relations for these two samples. We emphasize 
that the M91 and T04 abundance calibrations yield iden- 
tical conclusions. The red squares and blue circles show 
the median metalhcity in 0.05 and 0.15 dex wide bins of 
stellar mass for our SDSS and low-redshift AGES sample, 
respectively. The error bars on the AGES points refiect 
the standard error of the mean metalhcity in each stellar 
mass bin (the SDSS error bars are generally smaller than 
the symbols). 

Examining Figure 111! wc find that the two A4 — Z re- 
lations arc statistically consistent with one another. This 
figure also strikingly demonstrates that the SDSS sample 
is crucial for constraining the massive end of the Ai — Z 
relation; our low-redshift AGES sample simply contains 
too few massive, metal-rich galaxies due to the limited 
cosmological volume probed by the survey. Moreover, 
the consistency of the two Ai ~ Z relations suggests that 
the relative aperture bias between the two surveys is neg- 
ligible, even though the fibers used by AGES are a factor 
of two smaller (see also the detailed discussion of aper- 
ture effects in Therefore, we conclude that our 
SDSS sample provides an unbiased sample for measur- 
ing the evolution of the Ai — Z relation at intermediate 
redshift from AGES. 

6.1.2. B-band L — Z Relation at z ^ 0.1 

We turn our attention next to measuring the i3-band 
L — Z relation at z ~ 0.1 using our SDSS sample. Al- 
though the optical L — Z relation is considerably more 
sensitive to variations in star formation history (age or 
mass-to-light ratio, recent bursts of star formation, etc.) 
and dust attenuation than the Ai — Z relation (e.g., 
ILee et al.l[200l ISalzer et al.l[2005h . many previous stud- 
ies of the nebular abundances of star-forming galaxies at 
both low and high redshift have investigated the S-band 
L — Z relation as a proxy for the Ai — Z relation. How- 
ever, we will see in i 36.2.2l that the optical L — Z relation 
is an ineffective means of quantifying the chemical evolu- 
tion of star-forming galaxies due to its greater sensitivity 
to the effects of luminosity evolution. 

In Figure[T2]we show the i3-band L — Z relations from 
the SDSS based on the same four abundance calibrations 
presented in Figure [TO] As has been known for nearly 
three decades (see star- forming galaxies exhibit a 
strong correlation between optical luminosity and gas- 
phase abundance; however, the L ~ Z relation is clearly 
sensitive to the adopted abundance calibration, among 
other factors. To quantify the observed differences we 
defer to historical precedent and model the L — Z rela- 
tion using a simple linear model, given by 



12 + log (0/H) = Co + ci(Mb + 20.5), 



(11) 



where cq is the metalhcity at Mb = —20.5 and ci is 
the slope in dex mag~^. We solve equation pT|) using 
an ordinary least-squares bisector fit, which is appropri- 
ate when the scatter in the observed correlat ion is larger 
than the measurement error in either axis (jlsobe et al.l 
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Fig. 11. — M — Z relation at 2 ~ 0.1 based on our (left) SDSS and (right) low-redshift (0.05 < z < 0.20) AGES samples. In both panels 
the contours enclose 25%, 50%, 75%, 90%, and 97.5% of the sample, and the small grey points arc individual galaxies that lie outside the 
97.5% contour level. The open red squares and filled blue points represent the median mctallicity in 0.05 and 0.15 dex wide bins of stellar 
mass for the SDSS and AGES, respectively, and the error bars on the AGES points are the standard error in the mean mctallicity in each 
stellar mass bin. Where they overlap, the M — Z relations of the two samples are statistically consistent with one another. Note however, 
the relative paucity of massive, metal-rich galaxies in AGES over this redshift range, due to the relatively small amount of cosmological 
volume probed. 



ll99Clf ). We show the best-fitting L — Z relations in Fig- 
ure [T2] based on the KK04, M91, T04, and MPA-JHU 
metallicities as long-dashed red, dot-dashed blue, triple- 
dot-dashed green, and solid black lines, respectively, and 
list the best-fitting coefficients in Table |3l Again, the 
statistical uncertainties on the coefficients are negligible 
given the sample size and have not been tabulated. For 
comparison, we also show in the lower-right panel as a 
short-dashe d purple line the i3-ba. nd L — Z relation de- 
termined bv lTremonti et all (|2004l ) using the MPA-JHU 
abundances of ^ 53, 000 galaxies from the SDSS/DR4. 

We find that the slope of the S-band L — Z relation 
we derive ranges from —0.231 dex mag~^ when using the 
KK04 cahbration, to -0. 291 dex mag~^ based on the T04 
abundance calibration. iTremonti et al.l (j2004[ ). on the 
other hand, inferred a L — Z slope of —0.185 dex mag~^ 
using a sample that extended to much fainter luminosi- 
ties. Other studies studying both low-luminosity dwarfs 
and massive disk galaxies have obtained slopes rang- 
ing from —0.130 dex mag~^ to --0.280 dex mag~^ (e.g., 
Skillman et all Il989t iLee et al.l [20p: iLamareille et al.l 
2004 HSalzer' et al .l l2005tlMoustakas' et al We con- 

clude, therefore, that the slope of the L — Z relation is 
sensitive to both the adopted abundance calibration, and 
the range of luminosities spanned by the parent sample. 
Finally, we note that for a given abundance calibration, 
the la scatter in the U-band L — Z relation is roughly a 
factor of two larger than the scatter in the corresponding 
Ai — Z relation. 

6.1.3. Effect of Using EW(R2:i) on the Local M. — Z and 
L — Z Relations 



The preceding analysis has all been based on abun- 
dances derived from reddening-corrected fluxes; how- 
ever, to measure the evolution of the Ai — Z relation 
from AGES (see we must rely on oxygen abun- 

dances derived from emission- line EWs (see §5.2|1 . There- 
fore, in Table [3] we also list the coefficients of the 
SDSS M — Z and L — Z relations using our EW-based 
abundances. With respect to the Ai — Z relations de- 
rived using reddening-corrected fluxes, we find that the 
EW-based Ai — Z relations are slightly steeper (7 is 
smaller by 15% — 20%), but that the characteristic stel- 
lar masses and characteristic metallicities agree to within 
< 0.01 dex. These differences translate into maximum 
mctallicity differences of ^ 0.02 dex around ~ 10^'* A4q, 
and < 0.01 dex around 10^^ Mq. Similarly, the L — Z 
relations derived using either reddening-corrected fluxes 
or EWs have nearly identical slopes and zeropoints. 
Therefore, we conclude that using EWs to measure the 
evolution of the Ai — Z relation should not bias our con- 
clusions. 

6.2. Evolution of the M — Z and B-Band L — Z 
Relations Since z = 0.75 

Building on our measurement of the local Ai — Z and 
L — Z relations in the previous section, we now use 
AGES to measure their evolution since z = 0.75. 

6.2.1. Evolution of the M — Z Relation 

In Figure [T3] we plot the Ai — Z relation from AGES 
in six redshift bins centered on z = 0.1, 0.2, 0.3, 0.4, 0.5, 
and 0.65 (see also Fig. [T]). Here we only show results 
using the KK04 abundance calibration, although below 
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Fig. 12. — B-band L — Z relation from the SDSS based on four 
different abundance calibrations. The contours enclose 25%, 50%, 
75%, 90%, and 97.5% of the sample, and the small grey points rep- 
resent galaxies lying outside the 97.5% contour level. The dashed 
red, dot-dashed blue, triple-dot-dashed green, and solid black lines 
represent the ordinary least-squares bisector fits to the KK04, M91, 
T04, and MPA-JHU L — Z relations, respectively. We reproduce 
the MPA-JHU L — Z relation in each panel for reference, and in 
the lower-right panel we also show the L — Z relation derived by 
ITremonti et al.l l|2004l ) as a short-dashed purple line. The best- 
fitting slopes range from —0.231 dex mag~^ to —0.291 dex mag~^ 
(see Table [S]!, demonstrating that the choice of abundance cali- 
bration can have a significant effect on the infer red L — Z rela- 
tion. The large difference in slope with respect to ITremonti et aLl 
H2004l ^ . whose sample extended to lower absolute magnitudes, also 
demonstrates that the L — Z relation is sensitive to the range of 
luminosities spanned by the sample. 



we explicitly explore how the choice of abundance cal- 
ibration affects the amount of inferred evolution. For 
reference, we overplot the local M. — Z relation from 
the SDSS in every panel as a solid black curve (see Ta- 
ble [3]) . To quantify the observed evolution we compute 
the weighted mean metallicity in 0.2 dex wide bins of 
stellar mass, requiring a minimum of 10 galaxies with 
well-measured abundances in each bin. We weight each 
galaxy i by the effective inverse variance, given by 



_2 —2 

''cff,!: = ■^■»^log(o/H)i' 



(12) 



where Wi is that object's statistical weight (see ij2.2.ip . 
and o'iog(o/H)i is the statistical uncertainty on the metal- 
licity (see ij5.3p . We list the measured mean metallicities 
in Table m and plot the results in Figure [T^ using filled 
(open) red points for stellar masses above (below) our 
completeness limit (see i j4.3p . Examining Figure [iSl we 
find a well-defined M. ~ Z relation that changes remark- 
ably little over this rcdshift range. Nevertheless, at fixed 
stellar mass the mean metallicity of galaxies clearly de- 
creases with increasing redshift, although it is not yet 
clear whether the observed evolution varies with mass. 

We note in passing that the distribution of oxygen 
abundance at fixed stellar mass also appears to vary with 



redshift. For example, in the two highest redshift bins 
we see a tail of massive galaxies with lower-than-average 
oxygen abundances, 12-|-log (0/H) w 8.8—9, that are not 
present at lower redshift. We find a comparable number 
of objects based on both the T04 and M91 abundance 
calibrations. These objects may be similar t o the pop- 
ulatio n of "metal-poor" galax ies reported by iLillv et all 
()2003D and lMaier et al.l (|2005[) . However, we defer a more 
detailed analysis of this population of objects, and of the 
full metallicity distributions to a future paper, and focus 
instead on the evolution of the mean metallicity with 
redshift. 

To investigate the trends seen in Figure [T3] in more 
detail, in Figure [14] we show the mean metallicity (12 -|- 
log (0/H)) versus rcdshift for both our SDSS and AGES 
samples in eight bins of stellar mass spanning 10^'^ — 
10"-^ Mq. We show resuhs based on the KK04, T04, 
and M91 calibrations using red points, green triangles, 
and blue squares, respectively, and render stellar mass 
bins above and below our completeness limit using filled 
and open symbols. We model the observed evolution, 
separately for each calibration and stellar mass interval, 
as a linear function of redshift given by 



(12+log(0/H)) = (12+log(0/H)), 



d[log(0/H)] ^ 
=o.iH J [z~ 



(13) 

where (12 -I- log (0/H))z=o.i is the mean metallicity at 
z = 0.1 and d[log (0/H)]/dz is the logarithmic rate of 
metallicity evolution. Wc tabulate the best-fitting linear 
coefficients in Table [S] 

Figure [14] conveys several interesting results. First, 
we find that the KK04, T04, and M91 abundance cal- 
ibrations yield mean metallicities that differ in an ab- 
solute sense by significant amounts (no vertical offsets 
have been applied to the measurements). On average, 
the KK04-based metallicities are 0.093 dex (24%) and 
0.16 dex (45%) higher than the metallicities derived us- 
ing the T04 and M91 calibration, respectively. These 
systematic differences likely originate from the different 
assumptions and ingredients underlying the photoioniza- 
tion models used to cal i brate these strong-line methods 
fe.g.. iKewlev fc Dopital 120021 ). Understanding the ori- 
gin of these discrepancies is beyond the scope of this pa- 
per, although these results highlight that even within the 
class of theoretical abundance calibrations there are non- 
ne gligible systematic difi^e rences (see, e.g., the discussion 
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Despite the significant differences in the absolute abun- 
dances of the galaxies in our sample, all three calibra- 
tions reveal that star-forming galaxies of a given stellar 
mass become increasingly more metal poor with increas- 
ing redshift. The rate of metallicity evolution, however, 
clearly depends on the adopted abundance calibration 
(see Table [5]). Specifically, we find that d[log (0/H)]/dz 
based on the KK04 calibration is systematically shal- 
lower (by a factor of ^ 1.7 at > 3cr significance) than the 
slopes measured using the other two calibrations. We 
attribute this difference to the relatively fiat 7W — Z re- 
lation for massive galaxies implied by the KK04 calibra- 
tion (see Fig. [TU] top- left panel), which has the effect 
of reducing the rate of metallicity evolution we measure 
compared to the other two calibrations. From these data 
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TABLE 4 

Mean Oxygen Abundance in Bins of Stellar Mass and Redshift^ 



Redshift Range 






(log [M/M 


0))d (12 + log(O/H)KK04> 


(12 + log{0/H)M9i) 


(12 + log(O/H)T04) 


10.9 < log (M/Mq) < 11.1 


U.U5 




0.20 


1217 


0.14 


10.98 


9.095 ±0.010 


8.975 ±0.010 


9.055 ±0.010 


0.15 




0.25 


15 


0.20 


10.97 


9.099 ±0.010 


8.964 ±0.015 


9.057 ±0.020 


0.25 




0.35 


27 


0.32 


10.97 


9.055 ±0.018 


8.896 ±0.027 


8.976 ±0.029 


0.35 




0.45 


41 


0.39 


10.99 


9.067 ±0.010 


8.923 ±0.015 


8.997 ±0.019 


0.45 




0.55 


18 


0.50 


11.01 


9.039 ±0.010 


8.891 ±0.013 


8.965 ±0.016 


0.55 




0.75 


16 


0.62 


10.97 


9.032 ±0.021 


8.884 ±0.033 


8.952 ±0.031 


10.7 < log (M/Me) < 10.9 


0.05 




0.20 


4070 


0.12 


10.78 


9.102 ±0.010 


8.985 ±0.010 


9.073 ±0.010 


0.05 




0.15 


16 


0.12 


10.78 


9.106 ±0.010 


8.977 ±0.016 


9.072 ±0.019 


0.15 




0.25 


35 


0.22 


10.80 


9.074 ±0.010 


8.933 ±0.013 


9.015 ±0.015 


0.25 




0.35 


65 


0.31 


10.78 


9.077 ±0.010 


8.945 ±0.016 


9.025 ±0.018 


0.35 




0.45 


64 


0.39 


10.79 


9.065 ±0.010 


8.921 ±0.015 


9.004 ±0.015 


0.45 




0.55 


32 


0.51 


10.79 


9.039 ±0.010 


8.880 ±0.012 


8.963 ±0.014 


0.55 




0.75 


20 


0.61 


10.78 


9.058 ±0.015 


8.913 ±0.021 


8.968 ±0.040 


10.5 < log {M/Mq) < 10.7 


0.05 




0.20 


8760 


0.11 


10.59 


9.098 ±0.010 


8.980 ±0.010 


9.067 ±0.010 


0.05 




0.15 


19 


0.12 


10.58 


9.104 ±0.017 


8.982 ±0.022 


9.068 ±0.030 


0.15 




0.25 


75 


0.20 


10.58 


9.076 ±0.010 


8.938 ±0.017 


9.019 ±0.018 


0.25 




0.35 


84 


0.31 


10.60 


9.061 ±0.010 


8.910 ±0.013 


8.994 ±0.013 


0.35 




0.45 


80 


0.39 


10.59 


9.058 ±0.010 


8.909 ±0.011 


8.990 ±0.012 


0.45 




0.55 


31 


0.50 


10.61 


9.047 ±0.013 


8.898 ±0.020 


8.968 ±0.018 


0.55 




0.75 


14 


0.60 


10.63 


9.044 ±0.021 


8.885 ±0.031 


8.963 ±0.033 


10.3 < log (M/Me) < lU-5 


0.05 




0.20 


12679 


0.10 


10.40 


9.091 ±0.010 


8.963 ±0.010 


9.047 ±0.010 


0.05 




0.15 


37 


0.12 


10.39 


9.082 ±0.021 


8.963 ±0.041 


9.020 ±0.043 


0.15 




0.25 


120 


0.20 


10.39 


9.075 ±0.011 


8.937 ±0.020 


9.012 ±0.018 


0.25 




0.35 


114 


0.30 


10.39 


9.049 ±0.010 


8.899 ±0.012 


8.977 ±0.010 


0.35 




0.45 


93 


0.39 


10.39 


9.053 ±0.010 


8.899 ±0.011 


8.977 ±0.012 


0.45 




0.55 


28 


0.50 


10.39 


9.020 ±0.018 


8.847 ±0.026 


8.931 ±0.030 


10.1 < log {M/Mq) < 10.3 


0.05 




0.20 


14559 


0.08 


10.20 


9.076 ±0.010 


8.935 ±0.010 


9.014 ±0.010 


0.05 




0.15 


44 


0.11 


10.20 


9.072 ±0.010 


8.922 ±0.012 


9.007 ±0.014 


0.15 




0.25 


161 


0.21 


10.19 


9.053 ±0.010 


8.896 ±0.011 


8.974 ±0.012 


0.25 




0.35 


100 


0.30 


10.19 


9.037 ±0.010 


8.872 ±0.012 


8.952 ±0.012 


0.35 




0.45 


65 


0.39 


10.20 


9.025 ±0.010 


8.855 ±0.015 


8.925 ±0.016 


0.45 




0.55 


12 


0.50 


10.22 


8.965 ±0.014 


8.784 ±0.016 


8.854 ±0.019 


9.9 < iog(A^/A4Q) < 10.1 


0.05 




0.20 


13523 


0.08 


10.00 


9.053 ±0.010 


8.891 ±0.010 


8.967 ±0.010 


0.05 




0.15 


93 


0.12 


9.99 


9.032 ±0.011 


8.855 ±0.017 


8.924 ±0.016 


0.15 




0.25 


151 


0.21 


10.00 


9.017 ±0.012 


8.837 ±0.018 


8.903 ±0.021 


0.25 




0.35 


76 


0.30 


9.99 


9.012 ±0.010 


8.838 ±0.013 


8.919 ±0.013 


0.35 




0.45 


37 


0.38 


9.99 


8.953 ±0.015 


8.764 ±0.017 


8.832 ±0.018 


9.7 < log {M/Mq) < 9.9 


0.05 




0.20 


10378 


0.07 


9.81 


9.016 ±0.010 


8.833 ±0.010 


8.904 ±0.010 


0.05 




0.15 


105 


0.12 


9.80 


9.014 ±0.034 


8.827 ±0.035 


8.894 ±0.038 


0.15 




0.25 


162 


0.20 


9.80 


8.986 ±0.010 


8.797 ±0.012 


8.866 ±0.013 


0.25 




0.35 


63 


0.29 


9.81 


8.979 ±0.019 


8.780 ± 0.024 


8.850 ±0.025 


0.35 




0.45 


15 


0.37 


9.82 


8.955 ±0.022 


8.766 ±0.019 


8.833 ±0.018 


9.5 < log (M/Mq) < 9.7 


0.05 




0.20 


6474 


0.06 


9.61 


8.961 ±0.010 


8.763 ±0.010 


8.828 ±0.010 


0.05 




0.15 


107 


0.12 


9.60 


8.958 ±0.020 


8.764 ±0.022 


8.829 ±0.026 


0.15 




0.25 


118 


0.20 


9.59 


8.931 ±0.015 


8.736 ±0.014 


8.804 ±0.018 


0.25 




0.35 


33 


0.28 


9.62 


8.880 ± 0.020 


8.691 ±0.020 


8.747 ±0.023 



Mean metallicity of galaxies in multiple bins of stellar mass and redshift, based on the KK04, M91, and T04 calibrations. 
Note that the metallicities in the first row of each stellar mass interval correspond to our SDSS sample, while the metallicities 
in the other rows arc based on our AGES sample. 

^ Mean number of galaxies in this redshift interval and stellar mass bin. Note that the number of objects varies slightly 
depending on the abundance calibration (see, e.g., Table[TJ, so we show here the average number. We only provide redshift 
and stellar mass bins with at least ten galaxies. 

Mean redshift of the galaxies in this subsample. 

Mean stellar mass of the galaxies in this subsample. 
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Fig. 13. — AGES M — Z relation at z = 0.05 — 0.75 based on the KK04 abundance calibration. For reference, the contours enclose 25%, 
50%, 75%, and 90% of the galaxies in each redshift interval, and the small grey points represent individual galaxies that lie outside the 
90% contour level. The filled (opcn j po ints show the weighted mean metaUicity in 0.2 dex wide bins of stellar mass for galaxies above 
(below) our co mplete ness limit (sec q4.3|l . The solid black curve, reproduced in every panel for reference, is the local A4 — Z relation from 
the SDSS (sec ilG.l.lH . We find that galaxies experience relatively little chemical evolution over this redshift range. Note that the T04 and 
M91 abundance calibrations yield similar results. 
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TABLE 5 

Linear Evolution of the Mean Metallicity of Galaxies at z = 0.05 — 0.75^ 





(12 + log{0/H)>,=o.i 


d[log(0/H)]/dz 


(12 + log(0/H)),=o.i 


d[log(0/H)]/d2 


Calibration 


(dex) 


(dex z ^) 


(dex) 


(dex z ^) 






10.9 < yog{M/M 


0) < 11.1 


10.7 < log {M/M 


e) < 10.9 


KK04 


9.101 ± 0.01 


-0.139 ± 0.018 


9.104 ±0.01 


-0.140 ±0.015 


M91 


8.983 ±0.01 


-0.224 ±0.029 


8.990 ±0.01 


-0.236 ± 0.023 


T04 


9.064 ±0.01 


-0.244 ± 0.034 


9.079 ±0.01 


-0.276 ±0.028 




10.5 < [og[M/M 


a X 10.7 


10.3 < log [M/M 


q) < 10.5 


KK04 


9.100 ±0.01 


-0.161 ±0.021 


9.090 ±0.01 


-0.157 ±0.019 


M91 


8.982 ±0.01 


-0.284 ±0.033 


8.962 ±0.01 


-0.246 ±0.032 


T04 


9.069 ±0.01 


-0.306 ± 0.036 


9.045 ±0.01 


-0.279 ±0.031 




10.1 < [og[M/M 


0) < 10.3 


9.9 < log(A^/A^Q) < 10.1 


KK04 


9.073 ±0.01 


-0.177 ±0.031 


9.048 ±0.01 


-0.203 ± 0.040 


M91 


8.930 ±0.01 


-0.294 ±0.046 


8.884 ±0.01 


-0.267 ±0.054 


T04 


9.010 ±0.01 


-0.294 ±0.050 


8.960 ±0.01 


-0.256 ±0.056 




9.7 < log [M/M 


0j<9.9 


9.5 < log {M/M 


0)<9.7 


KK04 


9.010 ±0.01 


-0.225 ±0.076 


8.954 ±0.01 


-0.207 ±0.107 


M91 


8.825 ±0.01 


-0.267 ±0.088 


8.757 ±0.01 


-0.187 ±0.101 


T04 


8.895 ±0.01 


-0.286 ± 0.094 


8.822 ±0.01 


-0.159 ±0.128 



^ We model the measured change in the mean metallicity as a linear function of redshift given by: 
(12 + log (0/H)> = (12 + log (0/H)>,=o.i + d[log (0/H)]/d2 x (2 - 0.1), where (12 + log (0/H)),=o.i is 
the mean metallicity at z = 0.1 and d[log (0/H)]/d2: is the logarithmic rate of metallicity evolution. 
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Fig. 14. — Mean metallicity of galaxies at z = 0.05 — 0.75 as a 
function of stellar mass and redshift. The red points, green dia- 
monds, and blue squares correspond to the KK04, T04, and M91 
abundance calibration, respectively. Open and filled symbols rep- 
resent stellar masses above and below our completeness limit, re- 
spectively. We model the obs erve d evolution as a linear function 
of redshift given by equation I I13I I. and show the best fitting lines 
using solid red, dot-dashed green, and dashed blue lines for each 
calibration. We find that galaxies in all eight stellar mass bins 
become progressively more metal-poor with increasing redshift at 
comparable rates. 



alone it is impossible to ascertain which of these three 
strong-line calibrations is "right." The inevitable con- 
clusion, therefore, is that the uncertainties in both the 



absolute metallicites and the relative rates of metallicity 
evolution we measure are dominated by systematic errors 
in the strong-line abundance calibrations. 

With the preceding discussion in mind, we now ex- 
plore in more detail whether the chemical enrichment 
rate depends on stellar mass. In Figure [15] we plot 
d[log (0/H)]/dz versus stellar mass for galaxies with 
M = 10^-^ - 10" i Mq using the KK04 {red points), 
T04 (green triangles), and M91 (blue squares) abundance 
calibration. We exclude our lowest stellar mass bin in the 
following because it is most affected by the limited red- 
shift coverage of AGES in this stellar mass range. In each 
stellar mass interval we use the average mass of all the 
galaxies in that bin (i.e., across all redshifts), although 
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Fig. 15. — Redshift-dependent rate of change in the mean oxygen 
abundance for star-forming galaxies at 2 = 0.05 — 0.75 vs. stellar 
mass. The red points, green diamonds, and blue squares corre- 
spond to the KK04, T04, and M91 calibrations, respectively, while 
the corresponding solid red, dot-dashed green, and dashed blue 
lines show the best- fitting linear relation (see Table [Sj. The KK04 
calibration yields the smaller absolute rate of metallicity evolution, 
—0.160 dex z~^ at M = 10^"'^ A^0, but the steepest stellar mass 
dependence. By contrast, the T04 and M91 calibrations indicate a 
large rate of metallicity evolution, ~ —0.26 dex per unit redshift, 
and no significant dependence of this rate on stellar mass. 
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we have applied small horizontal offsets to the points in 
this figure for elarity. 

To quantify the observed trends we fit a linear model 
of the form 



d[log(0/H)] 
dz 



= ao + ai log 



M 



IQio-s Mr. 



(14) 



where ao is the rate of evolution in dex per unit red- 
shift, and ai eharaeterizes the power-law dependence of 
this rate on stellar mass. We list the best-fitting eoeffi- 
cients and uncertainties for the individual calibrations in 
Tabled 

Figure [15] nicely encapsulates many of the principal re- 
sults of this section, and of this paper. First, the factor of 
~ 1.7 shallower rate of metallieity evolution implied by 
the KK04 calibration compared with the T04 and M91 
calibrations is clearly apparent in this figure. Specifi- 
cally, the KK04 calibration indicates a metal enrichment 
rate of —0.160 ± 0.009 dex for star-forming galax- 
ies with M = 101°-^ 7^0, while the T04 and M91 cal- 
ibrations indicate more rapid evolution, —0.278 ± 0.015 
and -0.256 ± 0.014 dex respectively. The KK04 
calibration also indicates a steeper dependence of the 
rate of metallieity evolution on stellar mass, whereas 
the T04 and M91 calibrations show no statistically sig- 
nificant (< la) dependence on mass for galaxies with 
M. = 10^ ® — 10^^ Mq. We emphasize, however, that 
even the stellar mass dependence implied by the KK04 
calibration, 0.058 ±0.029, is only marginally statistically 
significant {2a). We conclude, therefore, that there is 
no compelling evidence that the M — Z relation evolves 
differentially, that is, in a mass-dependent way, over this 
range of stellar masses and redshifts. 

We conclude this section by examining whether we 
can synthesize the preceding results in terms of evolu- 
tion in the physical parameters 12 -I- log (0/H)* and Ai* 
given by equation ((TO)) (see also Appendix [B]). Recall 
that 12 + log (0/H)* is the asymptotic metallieity of the 
Ai — Z relation, and Ai* represents the stellar mass at 
which the relation begins to turn over. We begin with 
the simplest possible model by fixing Ai* and 7 at their 
local values and allowing 12 -I- log (0/H)* to vary linearly 
with redshift as 

12+log(0/H): = 12+log(0/H):^o.i+^(^-0.1), (15) 

where 12 + log (0/H)*^q is the characteristic metal- 
lieity at z = 0.1 (see Table and V is the rate of 
chemical evolution in dex per unit redshift. Solving 
for the maximum likelihood value of V using mpfitPH 
we obtain V = -0.137 ± 0.017, -0.258 ± 0.034, and 
-0.250 ± 0.044 dex z~i using the KK04, T04, and M91 
calibration, respectively. These results are statistically 
consistent with the rate of chemical evolution given in 
Table [HI (ignoring the weak dependence on stellar mass) , 
and therefore provide a convenient way of evaluating the 
A4 — Z relation for star-forming galaxies at any redshift 
between z = 0.05 and z = 0.75. 

In Figure [12] we show the evolution of the T04-based 
A4 — Z relation predicted by this empirical model, along 

A Lcvcnbcrg-Marquardt least-sq uares minimization rout ine 
available at http : //coh. physic s .wise . edu/$\slm$cralgm/idl[ 




Fig. 16. — Comparison of the T04-based mean metallicities of 
AGES galaxies in six redshift bin s betw een z = 0.1 and z = 0.6 to 
the empirical model proposed in i|6.2.1l The filled (open) symbols 
represent stellar mass bins above (below) our completeness limit in 
each redshift bin, and to best match the plotted measurements we 
have evaluated the model at the median redshift of the galaxies in 
each sample, as indicated in the legend. We find that this simple 
model, in which the J\4 — Z relation simply shifts toward lower 
metallieity without changing its shape, reproduces the observations 
well. 



with the mean metallicities from AGES in our six red- 
shift bins. As before, the filled and open symbols rep- 
resent stellar mass bins above and below our complete- 
ness limit, respectively. To optimally match the mea- 
surements plotted, we evaluate the model at the median 
redshift of the galaxies in each redshift bin, as indicated 
in the legend. We find that this simple model, which 
consists of a linear decrease in the oxygen abundance of 
star-forming galaxies as a function of redshift and an in- 
variant — Z relation shape, reproduces the observed 
data reasonably well. 

We also attempt to fit the data with a more complex 
model, in which the characteristic stellar mass Ml* also 
evolves linearly with redshift as 

log (MI/Mq) = log (Ml^o.i/Me) + Tl{z - 0.1). (16) 



Combining this relation with equation (|10|) , the redshift- 
dependent Ai — Z relation becomes 



12 + log (0/H), = 12 



log 



log (0/H) 1+7^(2 - 
M* X iO'R(-— "-1)^ 

109 TWq 



0.1) 

, (17) 



where once again we fix 7 at its local value. Finding 
the maximum likelihood values of the free parameters, 
we find that the data are only marginally better fit, but 
not by a statistically meaningful amount, compared to 
the simpler model in which 7?, = 0. We conclude that 
our metallieity measurements from the SDSS and AGES 
are statistically consistent with ?m Ai — Z relation that 
shifts toward lower metallieity with decreasing redshift, 
but whose shape does not evolve. 

6.2.2. Evolution of the B-Band L — Z Relation 
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Fig. 17. — Evolution of the B-band L — Z relation at z = 0.05 — 0.75 based on the KK04 abundance calibration. For reference, the 
contours enclose 25%, 50%, 75%, and 90% of the galaxies in each redshift interval, and the small grey points represent individu al gala xies 
that lie outside the 90% contour level. The solid black line reproduced in every panel shows th e loc al L — Z relation derived in ij6.1.2l and 
the vertical dashed lines show the completeness limits of AGES in each redshift interval (see i|4.3ll . Taken at face value, the observations 
suggest —1.00 dex of metallicity evolution at fixed luminosity {dashed red line), which significantly exceeds the much smaller rate of 
metallicity evolution inferred from the A4 — Z relation {dot-dashed blue line). We attribute the discrepancy to significant brightening in 
the _B-band luminosities of star-forming galaxies over this redshift range, and conclude that the i?-band L — Z relation does not provide 
an effective means of constraining the chemical enrichment histories of galaxies. 
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TABLE 6 

Stellar Mass Dependence of the Rate of 
Metallicity Evolution^ 



Calibration 


ao 






ai 




KK04 


-0.160 ± 


0.009 


0.058 


± 


0.029 


T04 


-0.278 ± 


0.015 


0.024 


± 


0.046 


M91 


-0.256 ± 


0.014 


0.052 


± 


0.042 



^ Mass-dependent rate of metallicity evo- 
lution given by d[log (0/H)]/d2: = ao + 
ai\og(M./lQ^^-^ -Mq)i where the units of ao 
are dex per unit rcdshift. 



We turn next to the evolution of the B-band L — Z 
relation at intermediate redshift. In Figure [IT] we plot 
the relationship between Mb and oxygen abundance in 
the same redshift bins used in the previous section. For 
clarity, we again only show the L — Z relation using the 
KK04 calibration, but we present results from the other 
two calibrations below. The solid black line reproduced 
in every panel is the SDSS L — Z relation derived in 
i j6.1.21 and the vertical dotted lines show the absolute 
magnitude limits of the survey in each redshift bin (see 
Table El). 

Taken at face value, Figure[T7]suggests that galaxies at 
z = 0.05 — 0.75 have experienced a significant amount of 
metallicity evolution, luminosity evolution, or both. To 
quantify the observed trends we modify the local L — Z 
relation given by equation to allow the characteristic 
metallicity at z = 0.1, co,z=o.ij to change linearly with 
redshift according to cq{z) = co,z=o.i + — 0.1), where 
the rate of metallicity evolution, S, is in units of dex z~^ . 
We also allow the ensemble of star-forming galaxies in 
AGES to brighten as Mb{z) = Mb,z=o.i + Qb{z - 0.1), 
where Qb is the rate of luminosity evolution in mag z~^. 
Combining these relations with equation (fTT]) . the B- 
band L — Z relation at a given redshift becomes 

12 + log (0/H)z = co,,=o.i + ci{Mb + 20.5) + 

(5 + ciQb)(2-0.1). (18) 

Note that we have implicitly assumed that the slope 
of the L — Z relation, ci , does not change with red- 
shift, which is likely an oversimplification; unfortunately, 
AGES does not extend far enough down the luminosity 
function to test this assumption. 

Momentarily ignoring luminosity evolution (i.e., set- 
ting Qb = 0), we find 5 -1.00, -1.32, and 
-1.11 dex using the KK04, T04, and M91 calibra- 
tion, respectively, with negligible statistical uncertain- 
ties. We show the results of this evolutionary model as 
a dashed red line in Figure [IT] for the KK04 calibra- 
tion. However, these results are clearly at odds with the 
~ —0.2 dex z^^ of evolution inferred from our analysis 
of the A4 — Z relation in j)6.2.11 which we plot in Fig- 
ure [17] dot-dashed blue line for comparison. We 
conclude, therefore, that the evolution of the i3-band 
L — Z relation requires star-forming galaxies in AGES 
to have been brighter in the past, that is Qb < 0. In- 
deed, setting V = S and solving for Qs, we obtain 
Qb = —3.5 ± 0.15 mag averaged over all three cali- 
brations. 

Measurements of the B-band luminosity function for 



blue, star-forming galaxies at intermediate redshift pro- 
vide an indep endent estimate of Qb- For example, 
iBlantonI ()2OO60 report a brightening of ^ 1.0 mag for 
blue galaxies since z = 1 using d ata fr om the SDSS 
and DEEP2 surveys; iFaber et al.1 ()2007f ) find a lumi- 
nosity evolution rate 1.4 mag z~^ over the same 
redshift interval from an analysis of COMBO- 17 and 
DEEP2 observations; and iCool et~all (|2011| ) find Qb = 
— 1.3 ± 0.2 mag z~^ based on data from the SDSS and 
AGES. Although we find a considerably larger rate of 
luminosity evolution, recall that the luminosity function 
studies cited above measure the brightening at the 'knee' 
of the luminosity function, M^, whereas our sample is 
dominated galaxies with Mb Mg aX higher redshift. 
Moreover, our assumption that the L — Z relation is lin- 
ear is clearly not strictly appropriate in light of our anal- 
ysis of the Ai — Z relation in §6.2.11 The principal con- 
clusion from the preceding analysis is that the i3-band 
L — Z relation does not provide a robust means of quan- 
tifying the evolution in the chemical abundances of star- 
forming galaxies at higher redshift, as ha s been previ- 
ously e mphasized by Tremonti et al.l ()2004[ ). ISalzer et aTl 
(pOOa ) , iZahid et all ([20111), and others. 



6.3. Comparison with Previous Results 

A detailed comparison of our results with previous 
studies poses numerous challenges. Among the issues 
to consider are the wide variations in sample selection, 
sample size, and survey area (i.e., cosmic variance), dif- 
ferences in the methods used to estimate oxygen abun- 
dances (i.e., choice of abundance diagnostics), and dif- 
ferences in the treatment of systematic effects like dust 
attenuation, stellar absorption, and AGN contamina- 
tion. Moreover, most previous abundance studies have 
measured the evolution of the i?-band L — Z relation, 
which we have seen is not a robust means of quantifying 
the chemical evolution histories of star-forming galax- 
ies. And finally, different studies have relied on disparate 
measurements of the local Ai — Z and L — Z relations, 
which can also infiuence the amount of chemical evolu- 
tion inferred. 

With these issues in mind, we focus our compari- 
son on three recent measureme nts of the Ai — Z rela- 
tion at intermed iate redshift by iCowie fc Barged (|2008l 
here after CB08) iLamareille et al.l (|2009l hereafter L09), 
and IZahid et all (|201lL hereafter ZKBll). We select 
these studies because they are based on a large num- 
ber of homogenously selected galaxies over a relatively 
wide area, and they make use of state-of-the-art tech- 
niques for measuring stellar mass from broadband pho- 
tometry. CB08 derived gas-phase oxygen abundances 
for roughly 200 galaxies at z = 0.05 — 0.9 in the 
145 arcmin^ GOODS-N (jGiavahsco et al.l [200l field, 
L09 measured the Ai — Z relation at z = 0.0 — 0.9 
based on optical spectroscopy of ^ 3000 galaxies in the 
0.61 deg2 VVDS-DEEP and 6.1 deg^ VVDS-WIDE sur- 
vey fields ()Le Fevre et al.ll2005t iGarilh et al.l[2008l ). and 
ZKBll measured the metallicities of 1350 galaxies at 
z = 0.75 - 0.82 selecte d from the 3.5 deg^ DEEP2 sur- 
vey (jDavis et al.l 120031 ) . For simplicity, we concentrate 
here on a single redshift interval, z = 0.6 — 0.8, for 
which all three studies have measured the Ai — Z re- 
lation. Although this redshift range is the highest red- 
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Fig. 18. — Comparison of the M — Z relation we infer at in- 
^rmediate redshift 
measurements by [Cowie 



termediate redshi ft from AGES ( f illed b lue squares) to previous 
Bargerf l l200g) ( ope n green diamonds) , 



ILamareille et al 



(open red circles), and lZahid et al.1 11201 ll ) 



(open pu rple triang l es), a ll converted to the T04 abundance scale 
and the IChabrierl 12003') IMF. For comparison, we show the 
M — Z relation at 2 ~ 0.1 (shaded black region) and at 2 ~ 0.6 
(shaded light blue regio n) based on the empirical evolutionary 
model derived in i|6.2.1l We find very good agreement between 
our empirical model, which is constrained at low stellar mass by 
our lower-redsh i ft obs ervations, a nd the detailed measur ements 
by IZahid et all ^OTH ). and with ICowie Barged ^OM ) above 
~ IQlO S Mq. 



shift bin probed by AGES, in which our sample spans a 
relatively limited range of stellar mass, the evolutionary 
model of the Ai — Z relation we derive in §6.2.11 is never- 
theless constrained by our larger sample of lower redshift 
galaxies. 

In Figure [THl we plot stellar mass versus oxygen abun- 
dance for galaxies at z ^ 0.7. We plot the mean 
metallicity measurements from CB08, L09, ZKBll, and 
AGES using green diamonds, red circles, purple trian- 
gles, and filled dark blue squares, respectively. For con- 
sistency across all the studies we show metallicities de- 
rived using the T04 calibration of R23. CB08 and L09 
derive their abundances using this calibration, and we 
have applied the approp riate conversion formula from 
iKewlev fc EUisoiil (|2008l ) to adjust the measurements 
from ZKBll. Finally, we subtract 0.25 dex from the 
CB08 stellar masses to convert them to the IChabrierl 
(|2003l) IMF. We also show in Figure HH the M-Z rela- 
tion at z ~ 0.05 — 0.15 as a shaded black region, and the 
M-Z relation at z = 0.50 - 0.70 from AGES as a hght 
blue shaded region. 

Examining Figure 1181 we find generally very good 
agreement between the intermediate-redshift M. — Z re- 
lation we infer from AGES, and the M — Z relations 
measured by CB08 and ZKBll. The agreement with 
ZKBll between - 10^-^ - 10^°-^ Mq is especially strik- 
ing since it represents an extrapolation of the A4 — Z 
relation we measure at lower redshift (AGES does not 
probe these stellar masses at this redshift). In partic- 
ular, our conclusion that the evolution of the Ai — Z 
relation does not depend significantly on stellar mass 
(see ij6.2.ip is strengthened by this comparison. The 
evolution of the Ai — Z relation since z ~ 0.7 is con- 



sistent with a progressive increase in metallicity, with 
no significant change in shape, contrar y to the con- 
clusions of several previous studies (e.g.. iSavaglio et ahl 
l2005t ILamareille et al.ll2"009( ). Our intermediate-redshift 
A4 — Z relation is also reasonably consistent with the 
measurements by CB08, particularly at the massive end. 
Between ~ 10^° - lO^^-^ Mq, however, CB08 measure 
slightly lower mean metallicities, but given the uncer- 
tainties in this comparison the differences are probably 
not significant. On the other hand, L09 measure mean 
abundances that are ^ 50% lower at ~ 10^" '^ Mq and 
a considerably shallower M — Z relation. The origin of 
this discrepancy is not clear, although the relatively low 
spectral resolution {R « 230) of the VVDS spectra ana- 
lyzed by L09 may introduce unintended selection effects 
related to a preferential loss of massive, metal-rich galax- 
ies with r elatively low equ ivalent-width emission lines 
(see, e.g.. iLiang et al.l[2004[ ). On balance, we conclude 
that our measurement of the M — Z relation at interme- 
diate redshift is consistent with previous measurements. 

7. SYSTEMATIC UNCERTAINTIES 

Before discussing our results in the next section, we 
first investigate how emission-line selection effects, aper- 
ture bias, AGN contamination, and other potential 
sources of systematic bias might impact our conclusions. 
Readers that are only interested in the interpretation of 
our results can safely skip ahead to ^ 

7.1. Emission-Line Selection Effects 

We selected our SDSS and AGES emission-line galaxy 
samples using an H/3 flux cut, and by requiring well- 
measured [O 11] A3727 and [O iii] A50G7 emission lines 
(see §4.ip . Here, we examine how these emission- line 
selection criteria might affect our conclusions. 

At a basic level, gas-phase abundances can only be 
measured for galaxies with ongoing star formation, oth- 
erwise H/3 will be absent from their integrated spec- 
tra. In principal, the metallicity of the gas during 
the last major epoch of star formation is encoded in 
the distribution of stellar metallicities of the stars re- 
sponsible for the continuum light; however, disentan- 
gling this relationshi p is complicated (alt hough not in- 
tractable; see, e.g.. ICaU azzi et al. 2005 : iToieiro et ahl 
l2007l:ICidFernandes et al.ll2007l:IPanter et al.ll2008f ). and 
outside the scope of the present analysis. At face value, 
therefore, our H/3 flux cut translates into an H/3 luminos- 
ity cut, or, modulo aper ture effects and du st attenuation, 
a cut in absolute SFR ()Kennicuttlll998b!) . 

In detail, however, our AGES and SDSS samples are 
also /- and r-band limited, respectively; therefore, H/3 in 
a galaxy at a given redshift must be sufficiently bright 
relative to the underlying continuum to be detected and 
measured given the design of the spectroscopic survey 
(instrumental setup, flux limit, exposure time, etc.). In 
other words, in a broadband flux-limited survey, an H/3 
flux cut translates into a EW-limitcd sample as a func- 
tion of redshift (see Fig.H]). This EW(H/3) limit is related 
to the minimum specifi c SFR (iKennicutt et al.l 119941 : 
iBrinchmann et al]|2004f ) a galaxy must exceed to enter 
our sample, although, the minimum (redshift-dependent) 
specific SFR depends on many details, including aperture 
effects, dust attenuation, the frequency of starbursts, and 
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the amount of luminosity and SFR evolution in galaxies. 
Nevertheless, comparison of our results against theoret- 
ical models should be aware of the systematic loss of 
galaxies with low specific SFRs with increasing redshift. 

In addition to excluding galaxies with low specific and 
absolute SFRs, our H/3 flux cut may also exclude ex- 
tremely dusty galaxies. In an integrated spectropho- 
tometric survey of star-forming galaxies in the nearby 
Universe that incl uded a significant numbe r of in frared- 
luminous galaxies. iMoustakas fc Kennicuttl (|2006l ) found 
that galaxies with significant Ha emission but with H/3 
too buried in the stellar continuum to be measurable 
constituted < 10 % of the sample with Lm > 10^^ Lq 
(so-called LIRGs; iSanders fc Mirabell [19961) . Given the 
rapid increase in the incidenc e of LIRGs to z = 1 
(|Le Floc'h et al.ll2005t iRodighie ro et al.ll2010( ). our sam- 
ple may be missing the dustiest galaxies at each redshift. 
On the other hand, we do not expect the metallicity prop- 
erties of the small fraction of very dusty galaxies to be 
significantly different with respect to the galaxies that 
satisfy our II/3 selection. Nevertheless, future studies of 
distant emission-line galaxies that are Ha- or even Paa- 
selected would mitigate the selection effects due to dust 
attenuation. 

Next, we also required our sample to have a well- 
measured [O ii] emission line. In ^4.11 we found that 
galaxies that failed this cut constituted ^ 6% of the 
H/3-selectcd sample and were frequently dusty starbursts. 
In a v olume-limited sample of SDSS galaxies, lYan et aP 
(|2006f l found that galaxies with significant Ha emission 
but without [O ii] emission constituted ~ 15% of their 
sample, the overwhelming majority of which had AGN- 
like forbidden emission-line ratios, and the remainder 
likely dusty star-forming galaxies. Therefore, we con- 
clude that requiring a well-measured [O ii] emission line 
has a negligible effect on our results. Indeed, the high 
detection rate of [O ii], due to its intrinsic strength 
and relative insensitivity to variations in metallicity and 
excitation across a wide range of physical conditions, 
is why this line remains an ap pealing SFR diagnostic 
at both low and h i gh redshift (IMoustakas et al. l 120061 : 
IGilbank et aD[2?nTil : iMostek et al.ll20lTi r 

Finally, we also required our sample of emission-line 
galaxies to have a significant [O ill] A5007 emission 
line. Unlike [O ll], the strength of [O III] decreases 
rapidly with increa sing metallicity and decreasing ex- 
citation (jKewlev fc Dopita 2002,); therefore, an [O iii]- 
limited sample could be biased against massive, metal- 
rich galaxies. To estimate the magnitude of this potential 
systematic bias, we appeal to the T04-based — ^ re- 
lation plotted in the lower- left panel of Figure [TOl Unlike 
the T04 calibration of R23 (see Appendix El), the MPA- 
JHU mctallicitics do not explicitly require significant 
[O III] emission. We find that the T04 and MPA-JHU 
A4 — Z relations agree very well, and only begin to devi- 
ate from one another around M sa 10^°'^ -Mq. Further- 
more, the magnitude of the deviation, in the sense that 
the T04-based Ai — Z relation at large masses might be 
underestimated, is always < 0.05 dex to 10^^'^ -Mq. We 
conclude, therefore, that requiring well-measured [O ill] 
emission line in our AGES and SDSS samples does not 
significantly bias our conclusions. 



7.2. Strong-Line Abundance Calibration 

We have shown that the choice of strong-line abun- 
dance calibration has a significant effect on the abso- 
lute normalization of the Ai ~ Z relation and, to second 
order, its shape (see Fig. [TU)) . The KK04 abundance 
calibration, for example, yields oxygen abundances that 
are, on average, 15% — 40% higher than abundances de- 
rived using the T04 and M91 calibrations, and results in 
a markedly shallower Ai — Z relation. At fixed stellar 
mass, these differences translate into a factor of 1.5 — 2 
systematic variation in the rate at which the mean metal- 
licity of star-forming galaxies changes with redshift (see 
gEOand Fig. [HI). 

These results highlight two important points. First, 
great care must exercised when examining chemical evo- 
lutionary trends with redshift based on oxygen abun- 
dances derived using different strong-line calibrations. 
For example, the abundances of distant galaxies derived 
using the KK04 calibration should not, in general, be 
compared blind l y to the M — Z relation published by 
iTremonti et al.l ()2004[ ). And second, even differential 
trends in metallicity can depend on the adopted abun- 
dance calibrat i on. A n ai ialogous result wa,s fou nd by 
iBresohn et all (|2009[ ) and IMoustakas et al.l (|2010[ ). who 
showed that the slope of the radial abundance gradient 
in nearby disk galaxies depends on the adopted strong- 
line calibration. Obviously both these effects have wide- 
ranging implications for using abund ance measurements 
to constrain theor etical models (e.g. JDutton et aLll2010l : 
IDave et al.ll20llB . 

Various techniques hav e been devised to deal with 
these issues. For example. iKewlev fc EUisonl (|2008| ) pro- 
vide polynomial expressions to transform abundances 
derived using different strong-line calibrations onto a 
common abundance scale. Other authors have used a 
combination of abundances derived using electron tem- 
perature measurements and photoionization models to 
attempt to generate a self-consistent suite of abun- 
dance diagnostics based on different s trong-line ratios 
(jNagao et al.ll2006l : lMaiolino et al.ll2008| ). affhough these 
methods are susceptible to systematic errors and in- 
consistencies in the photoionization models and abun- 
dance measurements based on electron temperature es- 
timates. However, these calibrations should not ap- 
plied beyond the range of parameter space for which 
they were deve loped (see especially the discussion in 
iStasihskal I2010D . Moreover, these techniques implic- 
itly ignore the origin of the systematic differences be- 
tween abundance derived using theoretical and empiri- 
cal strong-line methods (see ^Moustakas et al. 2010, and 
references therein), which may introduce metallicity- 
dependent biases. These issues are particularly rele- 
vant when applying strong-line calibrations based on the 
strength of the nitrogen emission line [N 11] A6584 to infer 
the oxygen abundances of galaxies, such as [N 11] /[O 11] , 
[N 11] /Ha, and ([O iii]/H/3)/([N 11] /Ha) (also known 
as the N202, N2, and 03N2 diagnostics, respectively; 
IKewlev fc Ellisonll2008l and references therein). Because 
nitrogen is both a primary and secondary nucleosyn- 
thetic product, the N/0 abundance ratio in galaxies 
has a well-known second-order dependence on the re- 
cent star formation history, which can introduce a non- 
negligible systematic bias in the inferred metallicities 
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Fig. 19. — Effect of AGN contamination on metallicities derived using the R23 parameter. The grey, blue, and orange contours, enclosing 
25% , 50%, 75%, 90%, and 97.5% of the galaxies in each sample, represent galaxies classified as SF, SF/AGN, and AGN, respectively (see 
^7.3l for the definition of these classes). {Left) Inspection of the T04-based M — Z relation clearly demonstrates that applying an _R23-based 
abundance calibration to galaxies hosting AGN results in metallicities that are spuriously low. (Right) Residuals from the M — Z relation 
vs. distance, D, from ([N Il]/Ha, [O lll]/H/3) = (—0.45,-0.5) dex in the BPT diagram. The solid black, dashed blue, and dot-dashed 
red curves show the median abundance residuals as a function of D for the SF, SF/AGN, and AGN class, respectively. Wc find that the 
metallicity bias for AGN is always > 0.2 dcx, and can be as high as ~ 0.6 dcx; meanwhile, the mctallicity bias for SF/AGN is always 
< 0.1 dex. 



(iContini et al .112003 iPilvugin et al.ll200llYin at alll2007l: 
iPerez-Montero fc Continill2009HThuan at al.ll20'Tol ). Ex- 
ploring thasa issues in mora quantitative detail is beyond 
the scope of this paper, suffice it to say that we have 
avoided most of these problems by using three indepen- 
dent abundance calibrations and applied them to our full 
sample of SDSS and AGES galaxies consistently. 

7.3. Residual AGN Contamination 

In H4.2l we identified AGN among the AGES emission- 
line eaJaxi^susinga combination of the BPT diagram, 
the lYan et al.l ()2011t ) empirical diagnostic diagram, and 
a variety of complementary multiwavelength diagnostics 
from the X-ray to the radio. We conservatively esti- 
mated, however, that up to ~ 50% of the emission- 
line galaxies in our z > 0.4 sample may host composite 
(weak) AGN. Of course, the true fraction may be higher 
or lower depending on a variety of different factors, such 
as evolution in the luminosity and optical color distri- 
bution of galaxies hosting AGN, dilution of the AGN 
spectrum due to aperture effects, and so forth. Never- 
theless, oxygen abundances derived from the R23 param- 
eter are particularly susceptible to AGN contamination 
because line-ratios involving high-ionization lines such as 
[O iii]/H/3 are prefere ntially elevated in t he narrow-line 
regions of AGN (e.g.. lOsterbrock fc FerTan d 2006). For 
galaxies on the upper R23 branch, an unidentified con- 
tribution from an AGN to the integrated emission-line 
spectrum will cause R23 to be overestimated^ and there- 
fore the metallicity to be systematically underestimated 
(see ij5.2l and Appendix R)). 



We use our SDSS sample to quantify the effect of 
unidentified AGN on our i?23-based oxygen abundances. 
In Figure [H (left) wc plot the T04-based M- Z re- 
lation for three sets of objects: (1) galaxies classi- 
fied as star-forming using the BPT diagram {grey con- 
tours, hereafter SF galaxies); (2) objects classified as 
AG N using th e BPT diagram, but as star-forming using 
the lYan et al.l ()2011[ ) diagnostic diagram {blue contours, 
hereafter SF/AGN); and (3 ) objects classifi ed as AGN 
using both the BPT and the lYan et al.l ()2011[ ) diagnostic 
diagram {orange contours, hereafter AGN). In this exam- 
ple we adopt the T04 calibration, but the same conclu- 
sions would hold using any other i?23-based abundance 
calibration. Figure [T2] {right) plots the residuals from 
the A4 ~ Z relation against the logarithmic distance, D, 
from log([N Ii]/Ha) = -0.45, log([0 lll]/H/3) = -0.5 
in the BPT diagram, at the base o f the sequence of 
star-fo rming galaxies (see Fig. 2] and iKauffmann et al] 
l2003c[) . For example, a galaxy with log([N ii]/Ha) = 
+0.2 and log([0 iii]/H^) ='+0.6 would have D = 
^(0.2 -f- 0.45)2 + (0.6 + 0.5)2 ^ 0.66 dex. In essence, 
powerful AGN have large values of D because they lie 
far from the tight sequence of star-forming emission-line 
galaxies, while normal star-forming galaxies have small 
D values. The solid black, dashed blue, and dot-dashed 
red curves show the median abundance residuals as a 
function of D for the SF, SF/AGN, and AGN, respec- 
tively. 

Figure [19] clearly reveals why AGN and, to a lesser 
extent SF/AGN, must be removed from any i?23-based 
abundance study: the inferred metallicities are spuri- 
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Fig. 20.— Fraction of /-band light subtended by the AGES I'.'S 
diameter fiber-aperture for galaxies in four bins of stellar mass 
vs. redshift. The large solid black squares correspond to the median 
light fraction (for all objects) in fixed 0.05- wide bins of redshift, and 
the error bars indicate the interquartile range of the light fraction in 
each bin. The median light-fraction is a weak function of redshift, 
increasing from ~ 13% at 2 ~ 0.1, to ~ 20% above z > 0.25. 

ously low because the emission lines are not due to pho- 
toionization from massive stars. Quantitatively, we find 
a very strong correlation between D and the metallicity 
bias. Among AGN, the metallicity can be (mistakenly) 
underestimated by as much as ^ 0.6 dex for objects with 
D ^ 1.2, with a typical bias of ^ 0.2 dcx for most objects. 
Meanwhile, the metallicity bias for SF/AGN is consider- 
ably more modest, ^0.1 dex (the separation betwen the 
solid black and dashed blue curves in the right panel) 
over the full range of D values. A 50 — 50 mixture of SF 
and SF/AGN, therefore, would exhibit a mean metal- 
licity bias of < 0.05 dex. We conclude, therefore, that 
unidentified SF/AGN in our sample of AGES galaxies, 
even those at z > 0.4, should not significantly bias our 
measurement of the evolution of the A4 — Z relation. 

7.4. Aperture Bias 

Most star-forming disk galaxies in the local Universe 
exhibit radial abundance gradients, in the sense that 
the metallicity in t heir central regions i s higher than 
in their outskirts (iZaritsky et all 119941 : Ivan Zee et al.l 
119981 : iMoustakas et al.ll2010( ). Therefore, oxygen abun- 
dances derived from fiber-optic spectra, which target 
the inner, metal-rich regions of the galaxy, may be 
biased relativ e to abundances derived from integrated 
spectroscopy (IZaritsky et al."19 95l:lKochanek et aLlbOOH: 
iKewlev et all 120051 : IMoustaka s et al.l 120061 ). If the 
amount of light subtended by the spectroscopic aperture 
is a strong function of redshift then so-called aperture 
bias (or aperture effects) could mimic genuine chemical 
evolution. 

In Figure [20] we explore this question by plotting the 
/-band light fraction versus redshift for galaxies in our 
abundance sample divided into four bins of stellar mass, 
log(A^/A^o): 9 - 9.5 (blue points), 9.5 — 10 {open tan 
squares), 10 — 10.5 {filled purple squares), and 10.5 — 11 
{green diamonds) . We define the light- fraction as the flux 
measured from our unconvolved /-band mosaics inside a 
1'.'5 diameter aperture centered on each galaxy, divided 



by the total galaxy flux. The median light-fraction for 
the whole sample is 0.19, with an interquartile range of 
0.14 — 0.24; however. Figure [20l reveals that the light- 
fraction is indeed a (weak) function of redshift. The large 
solid black squares correspond to the median light frac- 
tion in 0.05-widc bins of redshift, while the error bars 
indicate the interquartile range in each bin. The me- 
dian light-fraction increases from ~ 13% a.t z ^ 0.1, and 
then remains approximately constant at ^ 20% above 
z > 0.25. For reference, the physical diameter subtended 
by the AGES fiber-aperture is 5.9 kpc at z = 0.25, rang- 
ing from 2.8 — 10.7 kpc between z = 0.1 — 0.7. 

Taking these results at face value, the relative mctallic- 
ities we measure in AGES should be insensitive to aper- 
ture bias between z = 0.25 — 0.75. Below z ~ 0.25 we 
expect aperture effects to be more severe; however, at 
these redshifts we rely on the abundances inferred from 
the SDSS spectra, which were obtained through a 3" di- 
ameter fiber-optic aperture. The median i-band light- 
fraction for our SDSS sample is 0.24, with an interquar- 
tile range of 0.18— 0.31. The physical diameter subtended 
by the SDSS fiber is 5.5 kpc at z = 0.1, ranging from 
2.9 — 9.9 kpc between z = 0.05 — 0.2. In other words, by 
combining our SDSS sample at low redshift with AGES 
at higher redshift the relative metallicity evolution we 
measure should be insensitive to aperture effects. Of 
course, because of the existence of abundance gradients, 
the absolute metallicities we measure may be susceptible 
to aperture bias at all redshifts. 

We can place these results on a more quantitative foot- 
ing by performing some simple simulations. For con- 
venience, in the following we adopt the notation Z = 
12 -I- log (0/H). The inputs to our simulation are the (as- 
sumed linear) abundance gradient Z{r) = ZQ + {dZ/dr)r, 
where r is the galactocentric radius, Zq is the metallicity 
at r = 0, and dZ/dr is the slope; and the radial surface- 
brightness p rofile, I{r) (x exp[— (r/rg)"'^/"], where n is the 
ISersid (|1968D index, and tq is a sc ale factor that is related 
to the half-light radius, rso (see lGraham fc Drivei|[2005l 
for the relationship). The Sersic index indicates the con- 
centration of the galaxy, with increasing n corresponding 
to a mor e centrally concentrated surface- brightness pro- 
file fe.g.. lGrahai^l200ll : lPeng et al.|[200^ . For example, 
n = 1 implies a pure exponential disk, while n = 4 cor- 
responds to a classical [dc Vauco uleurd l|1948D r^/** light- 
profile. 

At a given radius the enclosed light-fraction, F{r), is 
given by 



F{r) 



I{r)rdr 



(19) 



and the surface-brightness weighted metallicity is given 

by 



{Z{r)) (X 



lo Z{r)I{r)rdr 
J^I{r)rdr 



(20) 



The integrated metallicity, (.^int)i that would be mea- 
sured from a spectroscopic aper ture that captures 
all the light of the galaxy (e.g., iJansen et al.l 120001 : 
IMoustakas fc Kennicuttir2006D follows from equation (|20|) 
by letting r — >■ (X). 
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Fig. 21. — Simulated effect of aperture bias on the metallicity 
inferred from spec tra t hat enclose a fraction of the integrated light 
of the galaxy (see i|7.4l for a complete description of the parameters 
of the simulation). For a given light-fraction, < F{r) < 1, the ab- 
scissa gives the difference betwreen the metallicity one would infer 
from a spectrum that encloses F{r) of the total light of the galaxy 
and the integrated metallicity. The metallicity bias is positive be- 
cause disk galaxies exhibit radial abundance gradients; therefore, 
their centers are generally more metal-rich than their outskirts. 
The contours and small grey points correspond to simulated galax- 
ies with a range of concentrations (Scrsic indices), sizes, and abun- 
dance gradient slopes. For reference, the contours enclose 50%, 
75%, and 95% of the simulated galaxies and the points are objects 
that lie outside the 95% contour level. The large blue squares show 
the median metallicity bias in fixed bins of F(r), and the error bars 
indicate the interquartile range. We find that the metallicity bias 
is typically < 0.1 dex, and a relatively weak function of enclosed 
light fraction. 



With these simple ingredients in-hand we ean now es- 
timate the relative effects of aperture bias on the in- 
ferred metallicities for a wide range of surface-brightness 
profiles, sizes, and abundance gradients. To constrain 
the range of n and r^o we use the distribution of phys- 
ica l properties for late-type ga laxies recently compiled 
bv iBlanton fc Moustakasl ()2009l . see especially §3.1 and 
Fig. 8 in that review) for SDSS galaxies at z < 0.05. 
Specifically, we adopt a normally distributed prior on n 
with mean (n) = 1.8 and standard deviation cr„ = 0.6, 
and a log-normal prior on r^o with (log {r^o)) = 0.35 kpc 
and criogf^rno) ~ '-'•25 kpc. For the range of abun- 
da nce gradient slope s we u se the recent measurements 
by iMoustakas et al.| ()2010f ) based on a sample of 21 
nearby disk galaxies; they find typical values (dZ/dr) = 
—0.03 dex kpc^^ and a^z/dr ~ 0-02 dex kpc~^. For 
each mock galaxy, we simulate the effects of metal- 
licity scatter at fixed galaetocentric radius (e.g., due 
to azimuthal metallicity variations) by perturbing the 
metallicity at each radius by a Gaussian distribution 
with mean 0.06 dex an d standard deviation 0.02 dex 
(|Moustakas et al.l |2010[ ). Finally, for each simulated 
galaxy we compute, as a function of r, the enclosed 
light fraction, F{r), and the metallicity bias, A{Z{r)) = 
{Z{r)) — (Zint), or the difference between the surface- 
brightness weighted metallicity at that radius and the 
integrated metallicity using equations ([T^ and (PO)) . We 
perform 1000 Monte Carlo simulations and plot the re- 
sults in Figure [21] 

Figure [21] reveals that, for our fiducial set of simula- 



tions, the metallicity bias is typically < 0.1 dex; more- 
over, it is a relatively weak function of enclosed light 
fraction. The principal reason for this result is that, 
although disk galaxies do indeed possess radial abun- 
dance gradients, they are not especially steep in abso- 
lute terms (i.e., in dex kpc~^). For example, making the 
mean abundance gradient slope a factor of three steeper, 
{dZ/dr) = —0.09 dex kpc~^ (e.g., if abundance gradi- 
ents were steeper at higher redshift; iBoissier fc Prantzoi 
l200l[ ). increases the absolute metallicity bias to 0.2 — 
0.3 dex for F{r) < 0.5, and makes it a much stronger 
of light-fraction; however, the relative metallicity bias 
over the relevant range of light-fractions (see Fig. [^0]) is 
still < 0.1 dex. Finally, we tested making our simulated 
galaxies a factor of 2 — 3 larger and a factor of 2 — 3 more 
concentrated and the basic results do not change. Based 
on these tests, we conclude that aperture bias should 
not significantly affect our measurement of the relative 
change in the abundances of star-forming galaxies from 
z = 0.05 - 0.75. 

8. DISCUSSION 

We have measured the evolution of the A4 — Z rela- 
tion from z ~ 0.05 — 0.75 using statistically complete 
samples of star- forming galaxies selected from AGES and 
the SDSS, carefully accounting for both random and sys- 
tematic sources of uncertainty. Our analysis has yielded 
three key results. First, we find that star-forming galax- 
ies obey a well-defined Ai — Z relation over the entire 
redshift range studied, z — 0.05 — 0.75. This result is 
generally consistent with previous measurements of the 
M. — Z relation at intermediate redshift, although the 
size and statistical completeness of our sample places it 
on a firmer quantitative footing. Second, we find that at 
fixed stellar mass, M ~ 10^°'^ -Mq, the mean metallic- 
ity of star-forming galaxies changes by just 30% — 60% 
since z ^ 0.7. The factor of two uncertainty on the 
amount of metallicity evolution is due to the systematic 
differences in the abundances inferred using the three 
strong-line calibrations investigated in this paper (KK04, 
T04, and M91). Finally, we find no statistically signifi- 
cant evidence that the A4 Z relation evolves in a mass- 
dependent way over this redshift range for galaxies with 
7W~lO9-«-lO"7W0. 

Despite considerable effort in recent years, a thorough 
theoretical explanation of the physical origin and evolu- 
tion of the M ~ Z relation remains elusive. However, 
models now provide quantitative estimates not only of 
the A4 — Z relation, but also auxiliary relations between 
stellar mass, SFR, and gas fraction, that can be utilized 
in concert to test the underlying model assumptions. Al- 
though we postpone a detailed comparison of our results 
on the evolution of the A4 — Z relation with various the- 
oretical models to a forthcoming paper, a brief discussion 
o f one class m o dels is suggestive of what can be done. 

[Dave et all (|2011a[ ) have proposed a simple an- 
alytic framework for interpreting the M — Z rela- 
tion based on insights gained analyzing a suite of 
cosmological hydrodynamie simulat ions with differ- 
ent pr e scriptions for galactic winds (iDave et al.l 120061, 
2011bt [Qppenheimer fc Pavel [2008t iFinlator fc Davel 
2008t lOppenheimer et al.ll2010( ). Their central hypoth- 
esis is that most galaxies are in a slowly evolving equi- 
librium with the intergalactic medium (IGM), balanced 
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between the competing effects of accretion, outflows, and 
star formation. They derive an expression for the M. — Z 
relation given by 



y 



(21) 



where 77 (the mass loading factor) is a proportionality 
constant relating the outflow-driven mass-loss rate to the 
instantaneous SFR, and az is the ratio of the metallic- 
ity of the gas being accreted from the IGM relative to 
the g as-phase metallicity of the galaxy (jFinlator fc Dav^ 
l2008f ). Their simulations show that a significant frac- 
tion of the metals previously ejected from galaxies in 
star-formation driven outflows ar e reaccreted in a mass- 
dependent way {wind recyclin g; lOppenheimer fc Pavel 
120081 lOppenheimer et alJ[20To[ ). or az > 0. The resuh 
is an — ^ relation whose shape is predominantly de- 
termined by the fraction of accreted gas that forms stars 
(modulated by 77), while its evolution is govered by az, 
the relative enrichment level of the infalling gas. Their 
favored mod el based on momentum-driven galactic winds 
(|Murrav et al. 200,53) grossly reproduces both the shape 
and evolution of the M — Z relatio n from z = ~ 3 
(with some notable discrepancies; see iDave et al.ll2011bl 
for details). 

An interesting consequence of this model is that de- 
viations from equilibrium at fixed stellar mass should 
correlate with metallicity, SFR, and gas fraction. For 
example, galaxies that have recently accreted significant 
amounts of cold gas from the IGM should have slightly 
lower metallicities (via dilution, assuming az ~ 0), and 
should be forming stars at slightly higher rates (due to 
the renewed gas supply) compared to galaxies of the 
same stellar mass still at equilibrium. These results pro- 
vide a compelling physical expla nation for the so-calle d 
fundamental metallicity relation (jMannucci et al.ll2010D . 
which posits that galaxies occupy a three-dimensional 
plane linking stellar mass, gas-phase metallicity, and 
SFR. Specifically, the fundamental metallicity relation 
reveals that at fixed stellar mass, galaxies with high SFRs 
tend to be metal-poor, while galaxies that lie above the 
median M — Z relation generally have lower SFRs. 

This type of comparison can, and will, be done with 
our SDSS and AGES observations. By incorporating ac- 
curate SFRs and approximate gas fractions, the latter of 
which we intend to esti mat e by inverting the Kennicutt- 
Schmidt ([K ennicutt 19983) law (see, e.g.. Appendix A in 
iPeeples fc Shankar,.201ll )7our analysis will provide pow- 
erful new insights into the physical origin and evolution 
of the M — Z relation, and, more generally, on our the- 
oretical understanding of galaxy formation. 

We conclude this section with a brief discussion of pre- 
vious and ongoing efforts to measure the A4 — Z rela- 
tion at even higher redshift, z > 1. Abundance studies 
of star-forming galaxies at these redshifts are especially 
challenging because they require near-infrared spec- 
troscopy, which remains difficult and time-consuming 
to obtain from the ground. Nevertheless, building on 
a handfu l of pioneering abu ndance studies of z 2 
galaxi es (iPettini et al.1 119981 l200ll iKobulnickv fc Kool 
l2000t IShaplev et all 120041: ISwinbank et al. I I2004D . mea- 
surements of the M — Z relation for star-forming 
galaxies now exist for galaxies at 1 < z < 1.6 



(IShaplev et all [2005t iMaier et al.|[200l iLiu et ai] [200l 
iPerez-Montero et al.|[2009HQuevrel et al.ll2009f) rLvman- 

break and near-infrared select ed g alaxies at z '--^ 2 
(|Erb et al.l 2006: Havashi et al.ll2009h . a nd star-forming^ 
galax i es at redshifts as high a s z ^ 3.5 ([M aiolino et al 



20081: iMannucci et aP 120091 : iLemoine-Busserolle et al 
loiaiRichard et al.ll201lD . 



It is instructive to highlight two important, outstand- 
ing issues confronting observational abundance studies of 
star- forming galaxies at z > 1. First, the availability of 
different emission lines depending on the redshift range 
of the sample under investigation and the spectral range 
of the observations means that metallicities in general 
cannot be derived using the same strong-line calibration 
at all redshifts. For example, the abundances of star- 
forming galaxies at z = 1 — 2 have been based, in most 
cases, on the N2 and 03N2 diagnostics fsee iff^ e.g. . 
IShaplev et al.ll2005HErb et al.ll2006l : IQuevrel et al.ll2009l 
but see iPerez-Montero et al.ll2009l for application of the 
so-called 02Ne3 diagnostic), while the R 23 diagnostic has 
been used for galaxies at z = 3 — 3.5 ([Mannucci et al.l 
120091 : iRichard et al.l[20Tll ) . Consequently, efforts to com- 
pare the M. — Z relation for star-forming galaxies at dif- 
ferent epochs is complicated by the systematic differences 
among strong-line abundance calibrations (see the dis- 
cussion in Second, observations suggest that the 
physical conditions (ionization parameter, electron den- 
sity or interstellar pressure, prevalence of shocks, etc.) in 
star- forming galaxies at z > 1 may be dramatically dif- 
ferent compared to local galaxies ( Shaplev et al.l 120041: 
Liu et al. 2008: iBrinchmann et all I2OO8I: iLehnert et atl 



20091: iHainline et al.ll2009i r It remains to be seen, there 



fore, whether using strong-line abundance diagnostics 
that have been calibrated against H 11 regions in nearby 
galaxies introduces systematic errors in the metallicity 
measurements of high-redshift galaxies. 

Addressing these issues requires concerted observa- 
tional effort on two fronts. First, high-quality observa- 
tions of H II regions in nearby galaxies, coupled with 
more detailed photoionization models, are needed to un- 
derstand the physical origin of the factor of ^ 5 sys- 
tematic uncertain ty in the nebular abun dance scale (see 
the discussion in iMoustakas et al.l 120101 and references 
therein). In addition, the effect of secondary nitrogen 
enrichment on oxygen abundances derived using the N2 
and 3N2 abundance d iagnostic s must be better under- 
stood (lYin et al.l [20071 : iPerez- Montero fc Conti^ 120091 : 
iThuan et al.ll2010i r And second, dedicated optical and 



near-infrared spectroscopic observing campaigns of large 
samples of higher redshift galaxies are needed to mea- 
sure the full suite of rest-frame optical emission-line di- 
agnostics. In particular, the advent of highly multiplexed 
near-infrared spectr ographs on 8 — 10-m class telescopes 
like LBT/lucifer (iMandel et all I2000D. Keck/MOSFiRE 
(IMcLean et"aIll2008D . GTC/emir (|Garzon et al.l l2003l) 
will facilitate the kinds of detailed studies that until re- 
cently have only been possible for lower-redshift galaxies. 

9. SUMMARY 

We have measured the gas-phase oxygen abundances, 
stellar masses, and i3-band luminosities of ^ 3000 star- 
forming galaxies at z = 0.05 — 0.75 observed as part 
of AGES, a redshift survey of Iab < 20.45 galaxies 
over 7.9 deg^ in the NDWFS Bootes field. This sam- 
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pie is among the largest statistically complete, wide-area 
samples of intermediate-redshift galaxies with measured 
nebular abundances assembled. We use state-of-the-art 
techniques coupled to high-resolution population synthe- 
sis models to measure the nebular emission lines free 
from the systematic effects of underlying stellar absorp- 
tion. In addition, we model the observed deep, aperture- 
matched broadband optical and near-infrared photome- 
try of the galaxies in our sample using iSEDf it, a new 
Bayesian SED-fitting code, to infer their stellar masses. 
Using multiple complementary multiwavelength diagnos- 
tics based on optical line-ratios and ancillary X-ray, mid- 
infrared, and radio observations, we identify and remove 
AGN from our sample. We combine volume-limited ob- 
servations from AGES with a complementary, statisti- 
cally complete sample of 75,000 star- forming galax- 
ies at z = 0.05 — 0.2 selected from the SDSS, which 
we analyze in the identical manner. We use the joint 
SDSS and AGES sample to measure the evolution of the 
M. — Z and B-band L — Z relations between z = 0.05 
and z = 0.75 using three independent strong-abundance 
calibrations (KK04, T04, and M91) of the metaUicity- 
sensitive i?23 parameter. 

We divide our principal findings into two categories, 
lessons learned and quantitative results. Among the 
lessons we have learned: 

1. EWs provide a powerful, reliable means of inferring 
oxygen abundances from nebular emission lines 
free from the systematic effects of dust attenua- 
tion. They also have the added advantage of being 
measureable in spectra that have not been flux- 
calibrated, although it is critical that the emission 
lines, especially H/3, be corrected for underlying 
stellar absorption. 

2. A thorough culling of AGN is necessary to ensure 
that metallicities derived using the i?23 parameter 
are not spuriously underestimated; the metallicity 
bias for galaxies with a subdominant (< 30%) con- 
tribution from an AGN is typically < 0.1 dex, but 
can be ^ 0.2 dex for powerful AGN. 

3. We find that aperture effects are not a significant 
source of systematic error in the gas-phase metallic- 
ities derived from fiber-optic spectroscopic surveys 
like the SDSS and AGES, owing to the relatively 
shallow radial abundance gradients in disk galaxies. 

4. The choice of strong-line calibration significantly 
affects not only the absolute normalization and 
shape of the M — Z relation, but also the rate of 
chemical evolution inferred, and its stellar mass 
dependence. Consequently, we caution strongly 
against comparing the metallicity measurements 
across redshifts based on different emission-line ra- 
tios (see especially the discussion in m.'2\i . 

5. We propose a new, physically motivated paramet- 
ric model for the M. — Z relation given by equa- 
tion ((TU)) . This parameterization has the same 
number of free parameters as the more commonly 
used polynomial model but has whose asymptotic 
behavior is much more stable. 



6. We confirm that stellar mass provides a much 
less ambiguous measurement of metallicity evolu- 
tion than optical luminosity, which itself evolves 
strongly. The i?-band L — Z relation in particular 
does not provide a reliable means of constraining 
the chemical evolution of star-forming galaxies. 

And our principal quantitative conclusions are: 

1. We find no statistically significant evidence for 
evolution in the shape of the M. ~ Z relation to 
z = 0.75 for M ~ 10^-^ - 10" Mq star-forming 
galaxies, contrary to previous findings based on 
smaller, incomplete samples. 

2. We find only modest evolution in the overall nor- 
malization of the M. — Z relation, suggesting a 
metallicity increase of between 30% — 60% for M ~ 
109-8 _ iQii Mq galaxies since z = 0.75. The un- 
certainty in our measurement is entirely due to the 
choice of strong-line calibration. 

Our results indicate that massive star-forming galaxies 
at intermediate redshift are chemically evolved, implying 
that they synthesized the bulk of their metals at higher 
redshift, z > 1. Given their stellar masses, this result is 
perhaps not too surprising. On the other hand, a chal- 
lenge for theoretical galaxy formation models that incor- 
porate chemical evolution will be to explain not only the 
relatively slow chemical enrichment rate for star-forming 
galaxies over this redshift range, but also the lack of evo- 
lution in the shape of the Ai ^ Z relation. In a forth- 
coming paper we will combine our metallicity and stellar 
mass estimates with SFRs and approximate gas fractions 
to directly test th e predictions of several of these theoret- 
ical models (e.g.. iBrooks et al.|[2007t iDave et al.ll2011al: 
IPeeples fc Shankadl201lD . 

On the observational side, joint optical and near- 
IR spectroscopy of intermediate-redshift star-forming 
galaxies would be especially useful. These observa- 
tions would provide the full suite of rest-frame opti- 
cal emission-line diagnostics, from [O ll] A3727 to Ha 
and [N ii] AA6548,6584, thereby enabling optical AGN 
to be identified more securely, and allowing gas-phase 
metallicities to be estimated using multiple independent 
strong-line abundance calibrations. Finally, abundance 
studies of statistically complete samples of intermediate- 
redshift galaxies spanning a wider range of st ellar masses, 
10^ — 10^^ Mq, should be a top priority fe.g.. lZahid et al.l 
120 lit ). Follow-up high-resolution optical and near-IR 
spectroscopy of lower-mass galaxies preselected from the 
latest gener ation of wide-area galaxy redsh ift surveys 
like G AMA (iDriver et al.llMll) zCOSMOS ()Lillv et al.l 
[2001 . and PRIMUS (ICoil et al.ll20Tll) would be one path 
forward. These abundance measurements, coupled with 
detailed stellar masses and SFRs, would provide valuable 
insight into the star formation and chemical enrichment 
histories of star-forming galaxies over a significant frac- 
tion of the age of the Universe. 

This paper has benefited from insightful sugges- 
tions, discussions, and encouragement from James Aird, 
Roberto Asset, Michael Blanton, Alison Coil, Aleks 
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APPENDIX 

A. ABUNDANCE CALIBRATIONS 

For completeness and the convenience of the reader, we describe and compare the three strong-line calibrations 
adopted in this paper. For the KK04 calibration, we have 

12 + log (O/H)KK04,iowor = 9.4 + 4.65a; - S.llx^ - log (q) (0.272 + 0.547x - 0M3x^), (Al) 

and 

12 + log (O/H)KK04,uppor = 9.72 - 0.777X - 0.951x^ ~ 0.0723;^ - 0.811x'' 

- log (g) (0.0737 - 0.0713x - O.Ulx^ -t- 0.0373x^ - O.OSSx"), (A2) 

for galaxies on the lower and upper branch, respectively, where x = log (i?23)- The ionization parameter q in cm s~^ 
is given by 



log (g) = 32.81 - 1.1537/2 + z(-3.396 - 0.025y + O.UUy^) 

X [4.603 - 0.3119y - 0.1632/^ + z(-0.48 + 0.0271y + 0m037y^)]~\ (A3) 

where z = 12 -I- log (0/H) and y = log (O32) characterizes the hardness of the ionizing radiation field. The ionization 
paramet er is also frequently written as a dim ensionless quantity U = q/c , where c = 2.99 x 10^"^ cm s~^ is the speed- 
of-light ()Shieldj|1990l:lKewiev &: Dopitall2002[ ). Note that equations (|Xl|) - (|X3)) must be solved iteratively for both the 
ionization parameter and the oxygen abundance; convergence is typically achieve d in a handful of i t eration s. 

For the M91 ca libration, we adopt the parameterizations proposed by iKobulnickv et all ()1999() . but see 
iKuzio de Naray et al.. (,2004) for alternative parameterizations based on trigonometric functions. For the lower branch, 
we have 

12 4- log (0/H)M9i,iowcr = 12 - 4.944 + 0.767a; + 0.602a;2 - y(0.29 + 0.332a; - 0.331a;2), (A4) 
and for the upper branch we have 

12 + log (0/H)M9i,uppcr = 12 - 2.939 - 0.2x - 0.237a;2 - 0M5x^ - 0.0283x'' 

- y(0.0047 - 0.0221a; - 0.102a;2 - 0.0817a;3 - 0.00717a;''). (A5) 
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Fig. A1. — Relationship between R23 and oxygen abundance for the three theoretical strong-line abundance calibrations adopted in this 
paper: T04 (solid red curve); KK04 [short-dashed blue and dotted orange curve, corresponding to log(C/) = —3.5 and —2.5, respectively); 
and M91 (dot-dashed green and long-dashed purple curve, corresponding to log(0^2) = —1.0 and 0.0, respectively). 

And finally the T04 calibration is only defined for galaxies on the upper i?23 branch: 



12 + log (O/H)T04.uppor = 9.185 - 0.313.T - Q.2Mx^ - OMlx^ (A6) 

In Figure lAll we compare the three separate strong-line calibrations for two values of the ionization parameters U 
and O32. 

B. PROPOSED NEW FUNCTIONAL FORM FOR THE M- Z RELATION 

Previous studies have parameterized the M ^ Z relation using a quadratic or third-order polynomial (e.g., 
iTremonti at al.l [20041 : IKewlev fc Ellisoiil [2008V A polynomial model is unsat isfactory, however, beca use it can lead 
to an — Z relation that decreases with increasing mass (see, e.g.. Fig. 2 of IKewlev fc E llison' ^008l). This artificial 
turn-over in the polynomial M. — Z relation is especially problematic when the model is extrapolated beyond the range 
of stellar masses for which it was calibrated. Therefore, in ^6.1.1l we fit the local — Z relation with a new functional 
form, which we reproduce here for the convenience of the reader: 



12 + log (O/H) = 12 + log (O/H)* - log 



1 



VlO^TWo 



(Bl) 



Unlike the polynomial model, this parameterization encapsulates our physical intuition that the M. — Z relation should 
vary monotonically with stellar mass. In our proposed model, M* refers to the stellar mass where the M. — Z relation 
begins to bend or flatten at a rate that is controlled by the power-law slope, 7. For example, an — Z relation 
with 7 <C 1 remains steep well above A^*, while an Ai — Z relation with 7 S> 1 is essentially flat for Al > Ai* . The 
parameter 12 + log (O/H)* is the asymptotic metallicity of the Ai ~ Z relation, i.e., the metallicity for galaxies with 
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The form and asymptotic limit of our model is motivated by the closed-box equations of chemical evolution 
(jTalbot fc Arnettlll971l : iBinnev k, Merrifieldlll998[ ). If the SFR declines exponentially with a characteristic timescale 
r, il}{t) (X exp(— t/r), then it is straightforward to show that the gas-phase metallicity at time i, Z{t), is given by 

Z{t) = -y In [1 - (1 - R) (1 - cxp-*/^)], (B2) 

where y is the nucleosynthetic yield, and R is th e return fraction, the fraction of mass returned by stars to the interstellar 
medium by supernovae and stellar winds ()Pagel,1997,lFl Note that in the abs ence of inflow and outflow, the return 
fraction simply equals the baryonic gas fraction. In the limit t 3> t equation (|B2p becomes Z = — yln(i?), i.e., the 
gas-phase metallicity becomes a constant equal to a fraction of the yield, which is analogous to our characteristic 
metallicity, 12 + log'(0/H)*. 

In addition to equation (|Bip and the traditional polynomial parameterization of the M — Z relation, we also inves- 
tigat ed a smooth double power-law model, wh ich is commonly used to fit the X-ray and infrared luminosity functions 
(e.g. lAird et al.|[2008l : iRuiopakarn et al.|[2"01ClD . and a broken double power-law. Unfortunately, the double power-law 
model exhibits the same general behavior as the polynomial model (i.e., it turns over rapidly at large stellar mass), 
while the broken power-law has a sharp, artificial break with no physical basis. 
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